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FORWARD 


The  algorithms  included  in  this  package  allow  the  user  to 
compute  first  order  standard,  percentage  and  logarithmic 
sensitivity  functions,  and  second  order  sensitivity  functions  of 
vehicle  transient  models  in  the  time  domain,  steady  state  models 
and  vehicle  models  in  the  frequency  domain. 

This  manual  assumes  that  users  already  have  a basic  knowledge  of 
FORTRAN  programming  and  some  experience  with  the  VAX/VMS  opera- 
ting system.  The  programs  included  in  this  package  have  undergone 
extensive  testing  and  have  proved  to  provide  useful  sensitivity 
results . 

The  methodology  utilized  in  the  development  of  this  software 
package,  as  well  as  descriptions  of  the  algorithms  and  sample 
sensitivity  results  of  vehicle  handling  models  are  presented  in 
the  F^nal  Technical  Report  - Part  II. 
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1.  PROGRAM  INSTALLATION 


The  programs  come  on  5 floppy  disks.  A summary  of  the  floppy 
disk  contents  are  shown  in  the  table  below: 

Summary  of  Disk  Contents 


Disk 

Number  File  Name 

Program  Description 

1 INSTALL.COM 

Com  file  used  to  install  software 

LIB.COM 

Com  file  used  to  compile  subroutines 

TRANS . FOR 

Main  program  for  sensitivity  of 
transient  vehicle  models  written  in  the 
time  domain 

VDATA . FOR 

Subroutine  used  to  read  parameter  and 
variable  information  from  a user 
generated  data  file 

INPUT 1. FOR 

Subroutine  used  for  interactive  input  of 
simulation  parameters 

DRVT . FOR 

Subroutine  used  to  fit  a polynomial 
through  various  points  and  calculates 
the  first  and  second  derivatives  at  a 
given  point 

FCN.FOR 

Function  subroutine  required  by  IMSL 
routine  DGEAR 

FCNJ . FOR 

Dummy  function  subroutine  required  by 
IMSL  routine  DGEAR 

GAUSS . FOR 

Subroutine  for  solution  to  simultaneous 
linear  equations  by  Gaussian  elimination 

SMD2 . FOR 

Sample  model  subroutine  of  spring-mass- 
damper  system  showing  model  structure 
required  for  primary  parameters  and 
primary  variables 

SMD3 . FOR 

Sample  model  subroutine  of  spring-mass- 
damper  system  showing  model  structure 
required  for  primary  parameters  and 
secondary  variables 

SMD4 . FOR 

Sample  model  subroutine  of  spring-mass- 
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damper  system  showing  model  structure 
required  for  secondary  parameters  and 
secondary  variables 

MODEL. FOR 

Sample  model  subroutine  of  three  degree 
of  freedom  non-linear  transient  vehicle 
model 

RO LOVER. FOR 

Sample  model  subroutine  based  on  STI 
rollover  model 

MODEL. DAT 

Parameter  data  file  for  three  degree  of 
freedom  vehicle  model 

ROLOVER . DAT 

Parameter  data  file  for  STI  based 

rollover  model 

SMD2 . DAT 

Parameter  data  file  for  spring-mass- 
damper  system 

SMD3 . DAT 

Parameter  data  file  for  spring-mass- 
damper  system  with  secondary  variables 

SMD4 . DAT 

Parameter  data  file  for  spring-mass- 
damper  system  with  secondary  variables 
and  secondary  parameters. 

STEDY. FOR 

Main  program  for  sensitivity  of  a steady 
state  vehicle  model 

INPUT2 . FOR 

Subroutine  used  for  interactive  input  of 
simulation  data 

SSFCN. FOR 

Function  subroutine  required  by  IMSL 
subroutine  ZSPOW 

SSMODEL. FOR 

Steady  state  model  of  three  degree  of 
freedom  passenger  vehicle 

SSMODEL.DAT 

Parameter  data  file  for  three  degree  of 
freedom  steady  state  vehicle  model 

FREQ DOM 3 . FOR 

Main  program  for  sensitivity  of  vehicle 
model  in  frequency  domain 

FVDATA . FOR 

Subroutine  used  to  read  parameter  and 
variable  information 

INPUT3 . FOR 

Subroutine  used  for  interactive  input  of 
simulation  data 
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TRANSFER. FOR 

Calculates  the  transfer  function  matrix 
of  a linear  mechanical  system 

FMODEL. FOR 

Three  degree  of  freedom  linear  frequency 
domain  model  of  passenger  vehicle 

FVDAT . DAT 

Parameter  data  file  for  frequency  domain 
vehicle  model 

4 PLOTSENS . BAS 

BASICA  plotting  program  for  transient 
time  domain  sensitivity  results 

5 PROCOMM . EXE 

PROCOMM  is  a public  domain  communica- 
tions software  package  used  for  file 
transfer  from  an  IBM  PC  to  any  mainframe 
computer  system. 

PRODOC24 . ARC 

PROCOMM  documentation  stored  in  archival 
form 

PROCOMM . PRM 

PROCOMM  file 

PROCOMM. XLT 

PROCOMM  file 

PROCOMM. KEY 

PROCOMM  file 

PROCOMM. DIR 

PROCOMM  file 

PROCOMM. HST 

PROCOMM  file 

1 . 1 Hardware  and  Software  Requirements 

To  install  and  run  the  sensitivity  algorithms  the  user  must  have 
access  to  the  following  hardware  and  software: 

1.  VAX  Mainframe  Computer  with  VMS  operating  system  and 

FORTRAN  77  compiler 

2.  IBM  PC  with  connection  to  VAX  mainframe 

3.  Hewlett  Packard  7475a  Plotter  (Optional) 

4.  Version  9.2  of  the  Double  Precision  IMSL  Library 
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1 . 2 Loading  Software  onto  VAX 


The  source  code  for  the  sensitivity  algorithms  is  written  onto 
floppy  disks  1-3 . The  source  code  must  be  uploaded  onto  a VAX 
mainframe . 

To  upload  the  programs  the  user  should  log  onto  the  VAX  mainframe 
from  their  PC.  This  can  be  done  by  inserting  the  PROCOMM  disk 
into  drive  a:  and  typing: 

A>PROCOMM 


This  will  connect  the  user  to  the  mainframe  provided  the  PC  is 
connected  directly  to  the  VAX.  If  the  PC  is  connected  by  a modem 
then  the  user  must  dial  up  the  VAX  and  connect  through  the  modem. 
If  PROCOMM  fails  to  connect  to  the  VAX  the  user  should  check  to 
see  if  the  communication  parameters  (Comm  port,  baud  rate, 
parity,  ect.)  are  properly  set.  This  can  be  done  in  PROCOMM  by 
pressing  Alt-p.  The  system  manager  should  be  able  to  advise  you 
on  the  proper  settings  if  required.  After  the  communication 
parameters  are  correct  the  user  should  then  logon  to  the  VAX. 

It  is  suggested  that  the  user  create  a separate  directory  for  the 
sensitivity  algorithms;  this  can  be  done  using  the  following  VAX 
VMS  command : 


$ CREATE/DIRECTORY  [.SENSITIV] 


The  user  should  then  switch  to  this  directory  using  the  following 
command : 


$ SET  DEF  [.SENSITIV] 


The  user  will  have  to  switch  to  this  directory  each  time  he  logs 
onto  the  VAX  and  wishes  to  use  the  sensitivity  algorithms. 

Once  the  user  has  changed  directories  he  should  copy  the  programs 
from  the  floppys  to  the  mainframe.  This  can  be  done  using  the 
PROCOMM  communication  software.  To  do  this  the  user  should  go 
through  the  following  steps. 


1.  Insert  Disk  1 into  Drive  b 

2 . Type  the  VMS  command  $ KERMIT 

3.  Type  the  command  Kermit-32>SERVER 

4.  Type  Alt  - k to  obtain  file  transfer  menu 

5.  Select  option  2 - Send. 

6.  Enter  Filespec  B:*.* 
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The  PC  should  start  sending  the  contents  of  the  floppy  to  the 
VAX.  The  user  then  insert  Disk  2 into  Drive  b and  repeat  steps  4 
thru  6,  and  repeat  the  same  for  Disk  3.  Once  the  contents  of 
Disks  1-3  have  been  sent  to  the  VAX  the  user  should  exit  KERMIT 
by  typing  Alt  - k and  selecting  option  3,  Finish. 

The  user  should  then  install  the  sensitivity  programs  by  typing 
the  following  command. 


$ @ INSTALL 


This  will  compile  the  source  codes  and  insert  them  into  three 
separate  object  libraries.  This  step  may  take  several  minutes. 
After  the  routines  have  been  compiled  they  must  be  linked  to  form 
executable  versions.  This  can  be  done  using  the  following 
commands . 


$ LINK  TRANS, TRNS/ LIB, IMS LI BD/ LIB 
$ LINK  STEDY, SDST/LIB, IMSLIBD/LIB 
$ LINK  FREQD0M3, FREQ/ LIB, IMSLIBD/LIB 


These  commands  will  create  the  following  executable  files: 


TRANS . EXE 
STEDY.EXE 
FREQD0M3 . EXE 


(Note:  If  problems  are  encountered  during  this  step  consult  your 
system  manager  about  the  command  syntax  required  to  link  to  the 
IMSL  library  on  your  system.) 

Once  this  step  has  been  completed  the  programs  are  in  executable 
form  and  ready  to  run. 
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2.  SENSITIVITY  ANALYSIS  IN  THE  TIME  DOMAIN 


****  TRANS. FOR  **** 


2 . 1 Program  Description 


Program  TRANS. FOR  calculates  first  order  standard,  percentage  and 
logarithmic  sensitivity  functions  as  well  as  second  order 
sensitivity  functions  of  transient  vehicle  models  written  in  the 
time  domain.  The  program  computes  the  parametric  sensitivity 
functions  defined  as  the  partial  derivative  of  system  variables 
with  respect  to  system  parameters.  The  system  must  be  modeled  as 
a system  of  differential  equations  (either  linear  or  non-linear) . 

The  program  stores  its  results  in  two  data  files.  The  first  file 
is  a user  readable  printer  output  data  file.  The  values  of  user 
selected  primary  and  secondary  variables  are  stored  in  this  file 
along  with  four  different  types  of  sensitivity  functions  and 
general  sensitivity  measures.  The  file  also  shows  the  order  of 
influence  of  system  parameters  on  system  variables  based  on  each 
type  of  sensitivity  function.  Hardcopy  of  this  file  Cun  be 
produced  by  sending  it  to  a line  printer.  The  second  data  file 
is  a plotter  output  data  file.  This  file  contains  variable 
values  as  well  as  the  four  different  types  of  sensitivity  func- 
tions and  general  sensitivity  measures.  The  plotter  output  data 
file  is  usually  stored  more  frequently  than  the  printer  output 
file  (for  better  graph  resolution)  and  is  stored  using  a file 
structure  which  is  more  compact  than  that  used  for  the  printer 
output  data  file. 

The  main  program  TRANS. FOR  uses  several  subroutines  including  the 
IMSL  routine  DGEAR,  there  is  also  a BASIC  program  used  to  plot 
output.  Short  descriptions  of  each  routine  are  given  below. 


DGEAR. OBJ  - A differential  equation  solver  which  is  part  of 
the  IMSL  library.  DGEAR  finds  numerical  solutions  to 
systems  of  first  order  differential  equations  with  given 
initial  conditions.  The  method  uses  a variable  step  size 
and  allows  the  user  to  select  a maximum  allowable  error 
criterion. 

FCN.FOR  - A function  subroutine  required  by  DGEAR. 

FCNJ . FOR  - A dummy  subroutine  required  by  DGEAR. 

INPUT1.F0R  - Provides  interactive  input  for  data  required  by 
the  sensitivity  algorithm.  This  information  includes  data 
required  by  the  IMSL  integration  routine  DGEAR,  the  time 
step  for  printer  and  plotter  output  storage,  the  output 
filenames,  the  total  simulation  duration,  and  initial  values 
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for  the  differential  equations. 

VDATA . FOR  - Reads  the  vehicle  parameters  from  a user 
generated  data  file.  This  file  includes  the  values  for  all 
primary  vehicle  parameters,  and  also  flags  which  variables 
and  parameters  (primary  and  secondary)  the  user  wishes  to 
include  in  sensitivity  calculations. 

MODEL. FOR  - Calculates  the  terms  which  are  integrated  by 
DGEAR  to  obtain  primary  variables,  also  calculates  secondary 
variables  and  secondary  parameters  needed  for  computing 
sensitivity  functions.  MODEL. FOR  is  interchangeable,  allow- 
ing the  sensitivity  algorithm  to  analyze  various  systems. 

DRVT . FOR  - Fits  a polynomial  through  any  given  number  of 
X,  Y coordinates  using  the  method  of  least  squares  and 
calculates  the  analytical  first  and  second  derivatives  of 
the  polynomial  at  the  desired  value  of  X. 

PLOTSENS . BAS  - A plotting  routine  written  in  BASIC  for  an 
IBM  PC  and  Hewlett  Packard  7475A  Plotter.  This  program 
reads  an  output  file  generated  by  the  sensitivity  algorithms 
and  allows  the  user  to  plot  graphs  of  variable  values, 
sensitivity  functions,  and  general  sensitivity  measures.  To 
utilize  this  program  the  user  must  to  transfer  data  from  the 
mainframe  computer  to  the  PC. 


2 . 2 Transient  Sensitivity  Model  Development 


In  order  to  analyze  the  sensitivity  of  a transient  time  domain 
model  the  user  must  develop  a model  according  to  certain 
guidelines  and  must  also  create  a parameter  and  variable  data 
file  using  a set  structure.  In  order  to  illustrate  the  steps 
involved  in  model  development  a sample  sensitivity  model  has  been 
developed  in  3 stages.  The  first  stage  analyzes  primary 
parameters  and  primary  variables  only.  The  second  stage  adds 
secondary  variables  and  the  third  stage  adds  secondary 
parameters.  It  is  strongly  suggested  that  users  developing 
models  for  the  first  time  implement  the  samples  given  below, 
before  attempting  to  develop  their  own  models. 

The  model  for  this  example  represents  a 2 degree  of  freedom 
spring-mass-damper  system  with  non-linear  damping  elements.  An 
illustration  of  this  model  is  shown  below: 
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Figure  1.  The  2 degree  of  freedom  spring-mass-damper  system. 


By  summing  the  forces  and  applying  Newtons  second  law  of  motion 
to  each  mass  the  following  two  simultaneous  second  order  dif- 
ferential equations  of  motion  are  obtained: 


m^Z]_  Z^(C^  t C2 ) — C2Z2  t + K2 ) Z2K2  — 0 


^2  Z 2 ( Z 2 “ Z ^ ) C2  ( Z 2 — Z-j  ) — F (t) 


1 1 ' 2 


1'  2 


The  model  is  then  brought  to  a first  order  system  of  differential 
equations  by  means  of  the  following  substitution: 

y1  = Z1 

y 2 = Z1 
y 3 = Z2 
y 4 = z2 
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This  substitution  reduces  the  system  to  a set  of  four  simul- 
taneous first  order  differential  equations: 


-1 

Yl  = [Yl(Ci  + C2)  ~ C2Y3  + Y2(k1  + K2 ) " Y4K2l 

ml 

Y2  = *1 

1 

Y3  = F(t)  - C(y3  “ Yl) c2  + ( Y 4 - Y2) k2] 

m2 

Y4  = y3 


After  the  transient  time  domain  model  has  been  brought  to  the 
first  order  form  it  is  possible  to  start  writing  the  model 
subroutine . 


2.2.1  Case  I - Primary  Parameters.  Primary  Variables 


The  simplest  form  of  the  model  subroutine  is  one  which  analyzes 
only  primary  parameters  (parameters  which  remain  constant)  and 
primary  variables  (variables  obtained  directly  from  integration) . 
The  first  example  shows  one  such  model.  The  vertical  coordinates 
Zi  = and  Z2  = Y4  of  masses  itu  and  m2  are  the  primary  variables 
of  interest,  and  the  spring  stiffnesses  K-^  and  K2  and  damping 
constants  ju  ^ and  ^ 2 are  the  parameters  of  interest. 

A print  out  of  the  required  model  subroutine  (File  SMD2 . FOR)  is 
shown  below: 


C *************************************************************** 


C * * 
C * SUBROUTINE  NAME:  MODEL  * 
C * * 
C * PURPOSE:  TO  CALCULATE  THE  RESPONSE  OF  A 2 DEGREE  OF  * 
C * FREEDOM  SPRING  MASS  DAMPER  SYSTEM  WITH  NON-  * 
C * LINEAR  DAMPING  COEFFICIENTS.  * 
C * * 
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* CALLING  SEQUENCE:  * 

* * 

* CALL  MODEL (N , X , Y , YPRIME , VAR, IFLG1)  * 

* * 

* ARGUMENTS : * 

* * 


* 

* 

* 

* 

* 

* 

* 

* 


N - NUMBER  OF  EQUATIONS  TO  INTEGRATE  * 

X - TIME  * 

Y - ARRAY  OF  STATE  VARIABLES  * 

YPRIME  - ARRAY  OF  INTEGRATION  ARGUMENTS  * 

VAR  - SECONDARY  VARIABLE  ARRAY  * 

IFLG1  - INTEGER  FLAG  (NOT  USED  IN  THIS  * 

EXAMPLE  BUT  SHOULD  BE  INCLUDED)  * 

* 


*•***************************************★*,********•******'****,*** 


C 

C 


c 

c 

c 


c 

c 

c 

c 

c 

c 


SUBROUTINE  MODEL (N , X , Y , YPRIME , VAR, IFLG1) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

REAL *8  Ml,  M2,  Kl,  K2 , MU1,  MU 2 , Cl,  C2 , Y(20),  YPRIME (20), 

! VAR ( 50 ) , Z10,  Z20 

COMMON  BLOCK  CONTAINING  PRIMARY  VEHICLE  PARAMETERS 

NOTE:  THE  ORDER  OF  PARAMETERS  IN  COMMON  BLOCK  VDAT  MUST  BE  IDENTICAL 
TO  THE  ORDER  IN  FILE  VDAT . DAT 

COMMON/VDAT/  Ml,  M2,  Kl , K2 , MU1,  MU 2 , Z10,  Z20 
CALL  VDATA 


CALCULATE  DAMPING  COEFFICIENTS 

Cl  = MU1  * DABS ( Y ( 1 ) ) 

C2  = MU2  * DABS ( Y(3)  - Y(l)  ) 


CALCULATE 

FORCE  OF 

SPRING  1 AND 

FORCE 

OF 

DASHPOT 

1 

FS1 

= ( Y(2) 

- Z10  ) * Kl 

FD1 

= Y(l)  * 

Cl 

CALCULATE 

FORCE  OF 

SPRING  2 AND 

FORCE 

OF 

DASHPOT 

2 

FS2 

= ( Y ( 4 ) 

- Z20  - Y ( 2 ) 

+ Z10 

) * 

K2 

FD2 

= ( Y ( 3 ) 

- Y ( 1)  ) * C2 

CALCLUATE 

EXCITATION  FORCE 

FI  = 0.0D0 

CALCULATE  EQUATIONS  OF  INTEGRATION 
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YPRIME ( 1 ) = -(  FS1  + FD1  - FS2  - FD2  ) / Ml 
YPRIME ( 2 ) = Y ( 1) 

YPRIME (3)  = FI  - ( FS2  + FD2  ) / M2 
YPRIME (4)  = Y ( 3 ) 

C 

C ADDITIONAL  CALL  TO  SUBROUTINE  VDATA 
C 

CALL  VDATA 

RETURN 

END 


The  transient  time  domain  sensitivity  algorithm  is  written  in 
double  precision  FORTRAN,  therefore  all  variables  must  be 
declared  as  double  precision.  The  current  version  of  the  algo- 
rithm is  limited  to  systems  represented  by  no  more  than  20  first 
order  differential  equations,  therefore  variables  Y and  YPRIME 
must  be  dimensioned  20.  The  algorithm  is  also  capable  of  exami- 
ning 50  different  variables  which  means  that  array  VAR  must  be 
dimensioned  50.  If  the  user  wishes  to  examine  systems  larger 
than  those  listed  above  then  he  must  increase  the  array  dimen- 
sions in  the  following  subroutines. 


TRANS . FOR 
FCN . FOR 
FCNJ . FOR 
MODEL. FOR 
INPUT 1. FOR 


The  common  block  VDAT  contains  all  primary  parameters.  It  is 
important  to  note  that  the  order  of  variables  in  common  block 
VDAT  must  be  identical  to  the  order  in  which  they  appear  in  data 
file  VDAT . DAT . More  information  will  be  given  regarding  the 
structure  of  data  file  VDAT.DAT  later. 

The  first  executable  statement  of  the  model  subroutine  must  be  a 
call  to  subroutine  VDATA.  This  allows  definition  of  the 
variables  in  common  block  VDATA.  After  this  call  has  been  made 
the  model  then  calculates  the  elements  of  array  YPRIME.  After 
the  elements  of  array  YPRIME  have  been  calculated  then  a second 
call  to  subroutine  VDATA  is  required.  After  this  has  been  done 
the  subroutine  is  finished  with  return  and  end  statements. 

The  file  VDAT.DAT  is  used  to  define  nominal  parameter  values, 
define  model  variables,  and  to  specify  which  parameters  and 
variables  the  user  wishes  to  use  in  the  sensitivity  analysis. 
Each  line  of  the  file  VDAT.DAT  must  be  written  in  a particular 
form.  This  form  is  shown  below: 
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0,'35  Characters  Used  for  Description  ',0.000 


The  first  field  of  the  line  is  an  integer  which  is  either  1 or  0 . 
The  second  field  is  35  character  long  and  is  used  for  descrip- 
tion. The  third  field  is  a real  number  of  double  precision  type. 
These  lines  are  added  together  to  form  the  mechanical  system's 
data  file.  One  possible  data  file  is  shown  below: 


0,  'Primary  Parameter  1 ',0.250D0 

1,  'Primary  Parameter  2 ',2.0D0 

1, 'Primary  Parameter  3 ',4.450D0 

0, ' ENDPP  ' , 0 . 0D0 

0, ' ENDSP  ' , 0 . 0D0 

0,  'Primary  Variable  1 Y(l)  ',0.0D0 

1,  'Primary  Variable  2 Y(2)  ',0.0D0 

0, 'Primary  Variable  3 Y(3)  ',0.0D0 

0, 'ENDOV  ' , 0 . 0D0 


The  first  lines  in  data  file  VDAT.DAT  are  used  to  denote  primary 
parameter  values.  This  particular  example  shows  three  primary 
parameters.  The  first  field  indicates  whether  sensitivity 
functions  will  be  calculated  with  respect  to  the  parameter.  A 
value  of  zero  signifies  that  sensitivity  functions  will  not  be 
calculated  with  respect  to  the  parameter  while  a value  of  one 
signifies  that  sensitivity  functions  will  be  calculated  with 
respect  to  the  parameter.  In  the  example  above  the  sensitivity 
functions  will  be  calculated  with  respect  to  the  second  and  third 
parameters  but  not  the  first.  The  second  field  represents  a 
description  of  the  parameter,  and  the  third  field  represents  the 
nominal  parameter  value.  As  mentioned  previously,  it  is 
important  that  the  parameters  are  listed  in  the  same  order  in 
which  they  appear  in  common  block  VDAT  in  the  model  subroutine. 

The  line  with  the  description  ENDPP  is  used  to  signal  the  program 
that  there  are  no  more  primary  parameters.  The  line  with  the 
description  ENDSP  is  used  to  signify  that  there  are  no  secondary 
parameters.  Even  though  there  are  no  secondary  parameters  in 
this  model  this  line  must  be  included  in  the  data  file. 

The  next  lines  in  the  data  file  are  used  to  describe  the  primary 
variables  of  the  model.  All  primary  variables  must  be  listed  in 
their  order  of  integration,  i.e.  Y(l),  Y(2),  ect.  The  first  field 
signals  whether  or  not  sensitivity  functions  of  this  variable 
will  be  calculated.  A value  of  0 signals  that  sensitivity  func- 
tions will  not  be  calculated  for  the  variable,  while  a value  of  1 
signals  that  sensitivity  functions  will  be  calculated  for  the 
variable,  therefore,  for  the  sample  file  above,  the  program  will 
calculate  the  sensitivity  functions  of  variable  2 with  respect  to 
parameters  2 and  3.  The  second  field  is  used  to  describe  the 
variable,  and  the  third  data  field  has  no  significance  and  should 
be  set  to  0.0D0. 
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The  last  line  in  the  file  has  the  description  ENDOV  which 
signifies  that  there  are  no  more  variables  and  that  the  end  of 
the  file  has  been  reached. 

The  sample  data  file  VDAT.DAT  for  the  spring-mass-damper  system 
described  previously  is  shown  below: 


0 , 'Mass  1 

0 ,  'Mass  2 

1 ,  'Spring  Stiffness  1 
1, 'Spring  Stiffness  2 
1, 'Damping  Factor  1 
1, 'Damping  Factor  2 

0, 'Relative  Position  Block  1 
0,' Relative  Position  Block  2 
0 , ' ENDPP 
0, ' ENDSP 


0, 

' 1 

Z1 

Dot 

= Y(l) 

If 

' 2 

Z1 

= Y(2) 

0, 

' 3 

Z2 

Dot 

= Y ( 3 ) 

If 

9 4 

Z2 

= Y ( 4 ) 

0 , ' ENDOV 


, 0. 5D0 
, 0.25D0 
, 10. 0D0 
, 5. 0D0 
, 1. 0D0 
, 0 . 50D0 
, 2 . 0D0 
, 4 . 0D0 
, 0 . 0D0 
, 0. 0D0 
, 0. 0D0 
, 0. 0D0 
, 0. 0D0 
, 0. 0D0 
, 0. 0D0 


Looking  at  the  data  file  we  can  see  that  the  program  will 
calculate  the  sensitivity  functions  of  variables  Y(2)  = and 

Y ( 4 ) = Z2,  which  represent  the  vertical  positions  of  masses  m-j 

and  m2  , with  respect  to  parameters  K2  , jm  ^ , and  p 2 . 

After  the  model  subroutine  has  been  written  and  the  file  VDAT.DAT 
has  been  generated  the  user  should  compile  the  model,  and  link 
the  routines  together.  The  model  SMD2 . FOR  should  be  compiled 

using  the  following  command: 


$ @LIB  TRNS  SMD2 


the  subroutines  are  linked  using  the  following  command: 


$ LINK  TRANS, TRNS/ LIB, IMS LI BD/ LIB 


Any  time  changes  are  made  to  a subroutine  the  subroutine  must  be 
recompiled  and  the  subroutines  linked  in  order  to  update  the 
executable  version  of  the  transient  sensitivity  algorithm.  Once 
this  has  been  accomplished  the  user  may  run  the  transient  sensi- 
tivity algorithm. 
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2.2.2 


Execution  of  Time  Domain  Sensitivity  Program 


Once  an  executable  version  of  the  program  is  obtained  and  the 
data  file  VDAT.DAT  is  created  the  program  can  be  executed  using 
the  following  command: 


$ RUN  TRANS 


The  program  then  gives  the  user  the  option  of  entering  the  data 
from  the  keyboard,  or  reading  it  from  an  existing  data  file.  The 
first  example  will  show  what  happens  when  the  data  is  entered 
from  the  keyboard. 


***  Data  Entry  *** 

1 . Enter  From  Keyboard 

2.  Read  From  File 

1 (For  input  from  keyboard) 

Input  Number  of  Equations  to  Integrate 
4 


Since  this  example  has  4 first  order  differential  equations  that 
must  be  simultaneously  integrated  the  user  should  enter  4 . 


Default  Values  For  Initial  Conditions 


Y( 

Y( 

Y( 

Y( 


1) 

2) 

3) 

4) 


0.0D0 
0 . 0D0 
0. 0D0 
0 . 0D0 


Do  you  wish  to  change  these  parameters  (Y/N) 
Y 


Since,  for  many  problems,  the  initial  values  of  the  differential 
equations  are  all  zero  the  program  uses  zero  as  the  default 
condition.  However  this  problem  has  non-zero  initial  conditions 
so  the  user  must  change  the  initial  conditions. 
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Input  Y ( 
0 . 0D0 

1) 

(Double 

Precision , 

i . e . 

0. 0D0) 

Input  Y ( 
0. 05D0 

2) 

(Double 

Precision, 

i . e . 

0.0D0) 

Input  Y ( 
0. 0D0 

3) 

(Double 

Precision , 

i . e . 

0 . 0D0) 

Input  Y ( 
0. 75D0 

4) 

(Double 

Precision, 

i . e . 

0. 0D0) 

Notice  that  the  initial  conditions  are  entered  as  double 
precision  real  numbers. 


***  Data  Required  By  IMSL  Routine  DGEAR  **** 

Default  Values 

H = . ID-05 
TOL  = . ID-09 
METH  =1 
MITER=2 
INDEX=1 


Do  you  wish  to  change  any  of  these  parameters  (Y/N) 
N 


These  numbers  represent  the  default  values  for  the  parameters 
required  by  the  IMSL  routine  DGEAR.  The  default  values  should 
prove  satisfactory  in  most  cases,  however  the  user  can  change 
these  values  should  it  prove  necessary.  The  IMSL  manual  should 
be  consulted  if  the  user  wishes  to  change  the  default  settings. 
For  this  example  the  default  values  will  not  be  changed. 


Input  the  Number  of  Sensitivity  Variables 
4 


This  represents  the  number  of  variables  which  can  be  included  in 
the  sensitivity  analysis.  This  includes  the  total  number  of  all 
primary  and  secondary  variables  listed  in  file  VDAT.DAT. 
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***  Program  Parameters  *** 

Derivative  Method  - Central  Difference  With  3 Points 
Percent  used  in  Central  Difference  - 1.0D-2 


Do  You  Wish  to  Change  These  Parameters  (Y/N) 
N 


Here  the  user  has  the  option  of  accepting  the  default  values 
specifying  the  number  of  points  used  for  derivative  evaluation 
and  the  parameter  change  associated  with  function  evaluations. 
The  user  may  select  a 3,  5,  or  7 point  central  difference  formula 
if  he  so  desires.  The  values  will  not  be  changed  for  this 
example . 


***  Simulation  Duration  and  Output  Increment  *** 

Input  Total  Run  Time  (Double  Precision) 

3 . 0D0 

Input  Printout  Time  Increment  (Double  Precision) 
0.5D0 

Input  Filename  For  Printer  Output  Data 
SMD2 . OUT 


Would  You  Like  the  Program  to  Create  a Plotting  Output  File  (Y/N) 

Y 


Input  Storage  Time  Increment  for  Plot  (Double  Precision) 
0.1D0 

Input  Filename  For  Plotting  Data 
SMD2 . PLT 


The  inputs  shown  above  indicate  that  the  sensitivity  functions 
will  be  calculated  for  the  first  three  seconds  of  the  model's 
response.  The  sensitivity  functions  will  be  written  to  the 
printer  output  data  file  every  0.50  seconds  and  the  printer 
output  file  name  will  be  SMD2.0UT.  A plotter  file  is  to  be 

created  and  will  be  stored  under  the  filename  SMD2.PLT.  The 

sensitivity  functions  will  be  stored  in  the  plotter  output  data 
file  every  0.10  seconds.  If  the  user  desires  to  create  a plotter 
output  file  then  the  time  increment  of  printer  output  should  be 
an  even  multiple  of  plotter  output  time  increment. 
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Do  You  Want  to  Store  This  Data  in  a File  (Y/N) 
Y 


What  Filename  Do  You  Wish  to  Give  the  Input  Data 
SMD2 . INI 


This  allows  the  user  to  store  the  information  entered  above 
(initial  conditions,  IMSL  data,  derivative  method,  simulation 
duration  and  data  output  information)  into  data  file  SMD.INI. 
Once  these  data  are  stored  in  a file  the  user  can  let  the  program 
read  the  data  from  the  file  on  subsequent  runs,  rather  than 
inputing  the  data  from  the  keyboard  each  time. 

After  all  of  the  above  questions  have  been  answered  the  program 
begins  calculating  sensitivity  functions  and  the  results  are 
stored  in  printer  and  plotter  output  files.  After  the  program 
has  stopped  the  user  can  send  the  results  stored  in  the  printer 
output  file  to  a line  printer  using  the  command: 


$ PRINT  SMD2.0UT 


Sample  results  obtained  after  the  first  printing  time  step  are 
shown  below: 


VARIABLE 

<N 

II 

> 

FUNCTION  VALUE  = 

0. 203356E+01 

1 

2 

TIME 

IP 

LDVDP-1 

DVDP-2 

0.50 

1 

-0. 138863E+00 

0. 342600E-02 

0.50 

2 

0. 144275E-01 

-0. 337260E-02 

0 . 50 

3 

0. 399735E-01 

-0.530420E-01 

0.50 

4 

-0. 283236E-02 

0. 105690E-01 

PARAMETER  EFFECT  ORDER  (1) : 1 3 2 4 

PARAMETER  EFFECT  ORDER  (2):  3 4 1 2 

SENSITIVITY  MEASURE  (1):  ALL  P ON  ONE  V=  0.14524783310 

SENSITIVITY  MEASURE  (2):  ALL  P ON  ONE  V=  5 . 5522444E-02 


VARIABLE  V=4  FUNCTION  VALUE  = 0.413169E+01 


TIME  IP 
0.50  1 


1 

LDVDP-1 
-0. 526390E-01 


2 

DVDP-2 

0. 119953E-02 
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0.50  2 

0.50  3 

0.50  4 


-0 . 725375E-01 
0.899959E-02 
0. 359551E-02 


0. 135132E-01 
-0. 184410E-01 
-0.256852E-01 


PARAMETER  EFFECT  ORDER  (1):  2 1 3 4 
PARAMETER  EFFECT  ORDER  (2) : 4 3 2 1 


SENSITIVITY  MEASURE  (1):  ALL  P ON  ONE  V=  9 . 014697121E-02 

SENSITIVITY  MEASURE  (2):  ALL  P ON  ONE  V=  3 . 440629534E-02 


The  VARIABLE  V=2  indicates  that  the  sensitivity  results  which 
follow  correspond  to  the  second  variable  listed  in  the  vehicle 
data  file  VDAT.DAT.  VARIABLE  V=4  indicates  that  the  sensitivity 
results  correspond  to  the  fourth  variable  listed  in  the  vehicle 
data  file.  The  value  of  each  variable  is  shown  next  to  the  label 
FUNCTION  VALUE. 

The  numbers  listed  under  the  column  labeled  IP  indicate  which 
parameters  correspond  to  the  sensitivity  values.  IP=1  represents 
the  first  parameter  selected  in  data  file  VDAT.DAT  (with  the 
first  data  field  set  equal  to  1)  , IP=2  the  second,  IP=3  the 
third,  ect. 

The  columns  labeled  LDVDP-l  and  DVDP-2  indicate  first  order 
logarithmic  and  second  order  sensitivity  values  of  variable  V 
with  respect  to  parameter  IP.  A total  of  four  different 
sensitivity  values  are  normally  stored  in  the  output  data  file, 
however  only  two  are  shown  in  this  example.  A description  of  the 
column  labels  are  shown  below: 


DV-l 

DVDP-1 

LDVDP-1 

DVDP-2 


First  Order  Percentage  Sensitivity 
First  Order  Standard  Sensitivity 
First  Order  Logarithmic  Sensitivity 
Second  Order  Standard  Sensitivity 


The  PARAMETER  EFFECT  ORDER  shows  the  order  of  influence  (from 
greatest  to  least)  of  all  parameters  on  variable  V based  on  each 
sensitivity  type.  The  SENSITIVITY  MEASURE  shows  the  combined 
effect  of  all  parameters  and  is  based  on  each  sensitivity  type. 
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2.2.3  Case  II  - Primary  Parameters.  Secondary  Variables 


Secondary  variables  are  variables  whose  values  are  not  obtained 
through  direct  integration  of  the  differential  equations  of 
motion.  Instead  they  represent  variables  which  are  functions  of 
system  parameters  and  primary  variables. 

If  analysis  of  secondary  variables  is  desired  then  model 
subroutine  and  system  data  file  must  be  modified. 

The  two  degree  of  freedom  spring-mass-damper  system  will  be  used 
as  an  example  to  demonstrate  the  modifications  required  by  the 
user.  The  forces  in  the  spring  and  damping  elements  supporting 
mass  m-^  will  be  used  as  secondary  variables  for  this  example.  A 
sample  listing  of  the  modified  model  subroutine  (File  SMD3.F0R) 
is  shown  below. 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


** 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

** 


★A***************************************************1******** 

* 


SUBROUTINE  NAME:  MODEL 

PURPOSE:  TO  CALCULATE  THE  RESPONSE  OF  A 2 DEGREE  OF 
FREEDOM  SPRING  MASS  DAMPER  SYSTEM  WITH  NON= 
LINEAR  DAMPING  COEFFICIENTS. 

CALLING  SEQUENCE: 

CALL  MODEL (N, X, Y, YPRIME, VAR, IFLG1) 

ARGUMENTS : 

N - NUMBER  OF  EQUATIONS  TO  INTEGRATE 

X - TIME 

Y - ARRAY  OF  STATE  VARIABLES 

YPRIME  - ARRAY  OF  INTEGRATION  ARGUMENTS 
VAR  - SECONDARY  VARIABLE  ARRAY 
IFLG1  - INTEGER  FLAG  (NOT  USED  IN  THIS 
EXAMPLE  BUT  SHOULD  BE  INCLUDED) 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 

* 

* 

★ 

* 

* 

* 

* 

* 


************************************************************* 


SUBROUTINE  MODEL (N , X , Y , YPRIME , VAR , IFLG1 ) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 


C 

REAL* 8 Ml,  M2,  Kl,  K2 , MU1,  MU2 , Cl,  C2 , Y(20),  YPRIME (20), 
! VAR ( 50 ) , Z10,  Z20 

COMMON  BLOCK  CONTAINING  PRIMARY  VEHICLE  PARAMETERS 
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c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


NOTE:  THE  ORDER  OF  PARAMETERS  IN  COMMON  BLOCK  VDAT  MUST  BE  IDENTICAL 
TO  THE  ORDER  IN  FILE  VDAT . DAT 

COMMON/VDAT/  Ml,  M2,  K1 , K2 , MU1,  MU2 , Z10,  Z20 
CALL  VDATA 


CALCULATE  DAMPINING  COEFFIECIENTS 

Cl  = MU1  * DABS ( Y ( 1)  ) 

C2  = MU2  * DABS ( Y(3)  - Y(l)  ) 


CALCULATE 

FORCE  OF 

SPRING  1 AND 

FORCE 

OF 

DASHPOT 

1 

FS1 

= ( Y(2) 

- Z10  ) * K1 

FD1 

= Y ( 1 ) * 

Cl 

CALCULATE 

FORCE  OF 

SPRING  2 AND 

FORCE 

OF 

DASHPOT 

2 

FS2 

= ( Y(4) 

- Z20  - Y ( 2 ) 

+ Z10 

) 

* K2 

FD2 

= ( Y ( 3 ) 

- Y ( 1 ) ) * C2 

CALC LU ATE 

EXCITATION  FORCE 

FI  = 0 „ 0D0 

CALCULATE  EQUATIONS  OF  INTEGRATION 

YPRIME ( 1 ) = -(  FS1  + FD1  - FS2  - FD2  ) / Ml 
Y PRIME (2)  = Y ( 1 ) 

YPRIME  (3)  = FI  - ( FS2  + FD2  ) / M2 
YPRIME (4)  = Y ( 3 ) 

DEFINE  SECONDARY  VARIABLES 

VAR ( 5 ) =FS1 
VAR ( 6 ) =FD1 

ADDITIONAL  CALL  TO  SUBROUTINE  VDATA 
CALL  VDATA 


RETURN 

END 


The  only  change  in  the  model,  as  compared  to  one  developed  prev- 
iously, is  the  addition  of  the  lines  defining  variables  VAR (5) 
and  VAR (6)  at  the  end  of  the  subroutine.  If  more  secondary 
variables  were  desired  then  they  could  be  defined  as  VAR (7) , 
VAR(8) , ect.  It  is  important  to  note  that  the  beginning  index 
for  array  VAR  is  5.  When  secondary  variables  are  defined  in  the 
model  subroutine  the  index  should  begin  at  a value  of  N+l  where  N 
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represents  the  number  of  primary  variables  (the  number  of  first 
order  differential  equations  of  motion.) 

Data  file  VDAT.DAT  must  also  be  altered  if  secondary  variables 
are  to  be  included.  A sample  of  VDAT.DAT  using  secondary  vari- 
ables is  shown  below: 


0, 

'Mass  1 

' , 0. 5D0 

0, 

'Mass  2 

' , 0. 25D0 

1, 

'Spring  Stiffness  1 

1 , 10.0D0 

1/ 

'Spring  Stiffness  2 

' , 5. 0D0 

1/ 

1 Damping  Factor  1 

' , 1. 0D0 

If 

' Damping  Factor  2 

• , 0. 50D0 

0, 

'Relative  Position  Block  1 

' , 2 . 0D0 

0, 

'Relative  Position  Block  2 

' , 4 . 0D0 

0, 

' ENDPP 

' , 0. 0D0 

0, 

' ENDSP 

' , 0. 0D0 

0, 

' 1 Z1  Dot  = Y (1) 

' , 0. 0D0 

If 

'2  Z1  = Y(2) 

' , 0. 0D0 

0, 

'3  Z2  Dot  = Y ( 3 ) 

' , 0. 0D0 

If 

'4  Z2  = Y (4) 

' , 0.0D0 

If 

' 5 Force  of  Spring  1 

' , 0. 0D0 

If 

' 6 Force  of  Dashpot  1 

' , 0. ODO 

0, 

' ENDOV 

' , 0. ODO 

The  only  change  in  this  file  is  the  addition  of  two  lines 
denoting  the  force  of  spring  1 and  the  force  of  dashpot  1 as 
variables . 

After  these  changes  have  been  made  to  the  model  subroutine  and 
data  file  the  user  can  then  compile  the  model,  link  the  program 
and  run  the  sensitivity  algorithm. 

The  model  should  be  complied  using  the  following  command: 


$ @LIB  TRNS  SMD3 


The  program  should  be  linked  using  the  following  command: 


$ LINK  TRANS , TRNS/ LIB, IMS LI BD/ LIB 


Program  execution  is  the  same  as  before  with  one  exception. 
After  inputing  the  initial  conditions  for  the  differential 
equations  of  motion  the  user  must  state  that  there  are  now  6 
sensitivity  variables  rather  than  the  previous  4 in  order  to 
account  for  the  two  additional  secondary  variables.  Output  ob- 
tained after  the  first  printing  time  step  is  shown-  below: 
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VARIABLE  V=2 

FUNCTION  VALUE  = 0.203356E+01 

TIME  IP 

0.50  1 

0.50  2 

0.50  3 

0.50  4 

1 

LDVDP-1 
-0. 138863E+00 
0. 144275E-01 
0. 399735E-01 
-0. 283236E-02 

2 

DVDP-2 

0. 342600E-02 
-0. 337260E-02 
-0. 530420E-01 
0. 105690E-01 

PARAMETER  EFFECT  ORDER  (1):  1 3 
PARAMETER  EFFECT  ORDER  (2):  3 4 

2 4 

1 2 

SENSITIVITY 

SENSITIVITY 

MEASURE  (1) : ALL  P ON 

MEASURE  ( 2 ) : ALL  P ON 

ONE  V=  0.14524783300 
ONE  V=  5.4297969E-02 

VARIABLE  V=4 

FUNCTION  VALUE  = 0.412169E+01 

TIME  IP 

0.50  1 

0.50  2 

0.50  3 

0.50  4 

1 

LDVDP-1 
-0. 526390E-01 
-0 . 725375E-01 
0. 899959E-02 
0. 359551E-02 

2 

DVDP-2 

0. 119953E-02 
0. 135132E-01 
-0. 184410E-01 
-0 . 256852E-01 

PARAMETER  EFFECT  ORDER  (1) : 2 1 
PARAMETER  EFFECT  ORDER  (2) : 4 3 

3 4 

2 1 

SENSITIVITY 

SENSITIVITY 

MEASURE  (1):  ALL  P ON 

MEASURE  ( 2 ) : ALL  P ON 

ONE  V=  9 . 014 69 7 IE- 02 
ONE  V=  3 . 4407055E-02 

VARIABLE  V=5 

FUNCTION  VALUE  = 0.335591E+00 

TIME  IP 

0.50  1 

0.50  2 

0.50  3 

0.50  4 

1 

LDVDP-1 
-0.7  4 14  06E+01 
0 . 874251E+00 
0 . 242225E+01 
-0. 171630E+00 

2 

DVDP-2 

-0.22217 IE-01 
-0. 337260E-01 
-0 . 530420E+00 
0. 105690E+00 

PARAMETER  EFFECT  ORDER  (1) : 1 3 2 4 
PARAMETER  EFFECT  ORDER  (2) : 3 4 2 1 


SENSITIVITY  MEASURE  (1): 
SENSITIVITY  MEASURE  (2): 


ALL  P ON  ONE  V=  7.850432538 
ALL  P ON  ONE  V=  0.542353017 
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VARIABLE  V=6  FUNCTION  VALUE  =-0 . 102447E+01 


1 

2 

TIME 

IP 

LDVDP-1 

DVDP-2 

0.50 

1 

-0. 262121E+00 

0. 160335E-01 

0.50 

2 

0.648263E+00 

-0. 162222E-01 

0.50 

3 

0 . 598784E+00 

0. 383520E+00 

0 . 50 

4 

0. 107629E+00 

0 . 284910E+00 

PARAMETER  EFFECT  ORDER  (1):  2 3 1 4 

PARAMETER  EFFECT  ORDER  (2) : 3 4 2 1 

SENSITIVITY  MEASURE  (1):  ALL  P ON  ONE  V=  0.926864889 

SENSITIVITY  MEASURE  (2):  ALL  P ON  ONE  V=  0.478311124 


2.2.4  Case  III  - Secondary  Parameters,  Secondary  Variables 


Secondary  parameters  are  parameters  which  do  not  remain  constant. 
Like  secondary  variables  they  are  functions  of  primary  parameters 
and  primary  variables,  in  fact  it  is  possible  to  have  a value 
which  is  treated  as  both  a secondary  parameter  and  secondary 
variable  simultaneously. 

The  two  degree  of  freedom  spring-mass-damper  system  will  again  be 
used  as  an  example  to  demonstrate  the  modifications  required  to 
investigate  the  influence  of  secondary  parameters  on  primary  and 
secondary  system  variables.  The  secondary  parameters  investigated 
in  this  example  are  the  damping  coefficients  of  both  dashpots  C-^ 
and  C2  , which  are  functions  of  the  primary  variables  Z-^  dot  and 
Z2  dot  and  primary  parameters  fUj.  and  ju2-  A sample  listing  of  the 
modified  model  subroutine  is  shown  below. 


C ********************** v**************************************** 


C * * 
C * SUBROUTINE  NAME:  MODEL  * 
C * * 
C * PURPOSE:  TO  CALCULATE  THE  RESPONSE  OF  A 2 DEGREE  OF  * 
C * FREEDOM  SPRING  MASS  DAMPER  SYSTEM  WITH  NON-  * 
C * LINEAR  DAMPING  COEFFICIENTS.  * 
C * * 
C * CALLING  SEQUENCE:  * 
C * * 
C * CALL  MODEL (N , X, Y , YPRIME , VAR, IFLG1 ) * 
C * * 
C * ARGUMENTS:  * 
C * * 
C * N - NUMBER  OF  EQUATIONS  TO  INTEGRATE  * 
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no  non  ooo  non  non 


c 

* 

X 

- TIME 

* 

c 

* 

Y 

- ARRAY  OF  STATE  VARIABLES 

* 

c 

* 

YPRIME 

- ARRAY  OF  INTEGRATION  ARGUMENTS 

* 

c 

* 

VAR 

- SECONDARY  VARIABLE  ARRAY 

* 

c 

* 

IFLG1 

- INTEGER  FLAG  (NOT  USED  IN  THIS 

* 

c 

* 

EXAMPLE  BUT  SHOULD  BE  INCLUDED) 

* 

c 

* 

* 

C *************************************************************** 


SUBROUTINE  MODEL (N , X , Y , YPRIME , VAR , IFLG1 ) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 


C 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


REAL* 8 Ml,  M2,  K1 , K2 , MU1,  MU 2 , Cl,  C2 , Y(20),  YPRIME(20), 

! VAR ( 50 ) , Z10,  Z20,  PCI,  PC 2 

COMMON  BLOCK  CONTAINING  PRIMARY  VEHICLE  PARAMETERS  ALONG  WITH 
PARAMETER  COEFFICIENTS  PCI  AND  PC2 

NOTE:  THE  ORDER  OF  PARAMETERS  IN  COMMON  BLOCK  VDAT  MUST  BE  IDENTICAL 
TO  THE  ORDER  IN  FILE  VDAT . DAT 

COMMON/VDAT/  Ml,  M2,  Kl,  K2 , MU1,  MU 2 , Z10,  Z20,  PCI,  PC2 

COMMON  BLOCK  CONTAINING  VALUES  OF  SECONDARY  PARAMETERS 

NOTE:  THE  ORDER  OF  PARAMETERS  IN  COMMON  BLOCK  SECP  MUST  BE  IDENTICAL 
TO  THE  ORDER  OF  PARAMETER  COEFFICIENTS  IN  FILE  VDAT . DAT 

COMMON/SECP/  Cl,  C2 

CALL  VDATA 


CALCULATE  DAMPINING  COEFFIECIENTS 


Cl  = MU1  * DABS ( Y ( 1 ) ) 

C2  = MU 2 * DABS ( Y(3)  - Y(l)  ) 

MULTIPLY  SECONDARY  PARAMETERS  BY  PARAMETER  COEFFICIENTS 

Cl  = Cl  * PCI 
C2  = C2  * PC2 


CALCULATE 

FORCE  OF 

SPRING  1 AND 

FORCE 

OF 

DASHPOT 

1 

FS1 

= ( Y(2) 

- Z10  ) * Kl 

FD1 

= Y(l)  * 

Cl 

CALCULATE 

FORCE  OF 

SPRING  2 AND 

FORCE 

OF 

DASHPOT 

2 

FS2 

= ( Y ( 4 ) 

- Z20  - Y (2) 

+ Z10 

) * 

K2 

FD2 

= ( Y ( 3 ) 

- Y ( 1 ) ) * C2 

CALCLUATE 

EXCITATION  FORCE 
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c 

FI  = O.ODO 
C 

C CALCULATE  EQUATIONS  OF  INTEGRATION 
C 

YPRIME ( 1)  = -(  FS1  + FD1  - FS2  - FD2  ) / Ml 
YPRIME (2)  = Y ( 1 ) 

YPRIME (3)  = FI  - ( FS2  + FD2  ) / M2 
YPRIME (4)  = Y ( 3 ) 

C 

C DEFINE  SECONDARY  VARIABLES 
C 

VAR ( 5 ) =FS1 
VAR ( 6 ) =FD1 
C 

C ADDITIONAL  CALL  TO  SUBROUTINE  VDATA 
C 

CALL  VDATA 

RETURN 

END 


There  are  substantial  differences  between  this  model  and  the 
models  shown  previously.  The  first  difference  is  the  addition  of 
two  variable  parameters  PCI  and  PC2  to  the  common  block  VDAT. 
These  two  variable  parameters  represent  secondary  parameter  coef- 
ficients. They  have  a nominal  value  of  1.0  and  are  varied  by  the 
percentage  of  parameter  variation  specified  by  the  user  to  obtain 
differing  model  responses.  These  secondary  parameter  coef- 
ficients PCI  and  PC2  must  be  multiplied  by  the  secondary  parame- 
ters of  interest  to  the  user.  The  secondary  parameter  coef- 
ficients PCI  and  PC2  must  be  included  in  the  data  file  VDAT.DAT 
after  the  ENDPP  marker  and  before  the  ENDSP  marker  as  shown 
below. 


0, 

0, 

1/ 

1, 

1, 

1, 

0, 

0, 

0, 

1, 

1, 

0, 

0, 

1/ 

0, 


Mass  1 

1 , 0. 5D0 

Mass  2 

! , 0. 25D0 

Spring  Stiffness  1 

' , 10. 0D0 

Spring  Stiffness  2 

* , 5. 0D0 

Damping  Factor  1 

' , 1. 0D0 

Damping  Factor  2 

' , 0. 50D0 

Relative  Position  Block  1 

' , 2 . 0D0 

Relative  Position  Block  2 

' , 4 . 0D0 

ENDPP 

' , 0. 0D0 

Secondary  Parameter  Coefficient 

PCI 8 , 1. 0D0 

Secondary  Parameter  Coefficient 

PC2 ' , 1. 0D0 

ENDSP 

8 , 0. 0D0 

1 Z1  Dot  = Y ( 1 ) 

' , O.ODO 

2 Z 1 = Y ( 2 ) 

' , O.ODO 

3 Z2  Dot  = Y ( 3 ) 

' , 0. 0D0 
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1 , ' 4 Z2 
1 , ' 5 Force 
1 , ' 6 Force 
0 , ' ENDOV 


= Y(4) 
of  Spring  1 
of  Dashpot  1 


' , 0. 0D0 
' , 0. 0D0 
' , 0. 0D0 
' , 0. 0D0 


The  next  difference  in  the  model  is  the  addition  of  the  common 
block  SECP . This  common  block  contains  the  values  of  the  second- 
ary parameters  selected  by  the  user  (Cl  and  C2)  and  is  used  to 
pass  the  computed  values  of  secondary  parameters  back  to  the 
sensitivity  algorithm.  These  secondary  parameters  must  be  placed 
in  the  common  block  having  the  same  order  as  the  corresponding 
secondary  parameter  coefficients  PCI  and  PC2  have  in  common  block 
VDAT. 

After  the  secondary  parameters  of  interest  have  been  written  in 
the  model  subroutine  the  secondary  parameters  should  be  mult- 
iplied by  their  corresponding  secondary  parameter  coefficients. 

Once  the  model  and  data  file  have  been  modified  to  accept 
secondary  parameters  the  model  subroutine  must  be  compiled  and 
linked  as  before.  The  procedure  for  running  the  program  is  the 
same  as  that  used  in  the  example  demonstrating  the  use  of  primary 
parameters  and  secondary  variables.  The  example  output  obtained 
after  the  first  printer  time  step  is  shown  below: 


VARIABLE 

CM 

II 

> 

FUNCTION  VALUE  = 

0. 157174E+01 

1 

2 

TIME 

IP 

LDVDP-1 

DVDP-2 

0 . 50 

1 

0 . 451245E+Q0 

-0. 729950E-02 

0.50 

2 

-0 . 282 153E-01 

0. 150981E-02 

0.50 

3 

-0. 179749E+00 

0. 258233E+00 

0 . 50 

4 

0. 237402E-01 

-0. 176348E+00 

0 . 50 

5 

-0 . 274720E+00 

0. 441750E-01 

0.50 

6 

0. 346095E-01 

-0.314228E-01 

PARAMETER  EFFECT  ORDER  (1) : 1 5 3 6 2 4 

PARAMETER  EFFECT  ORDER  (2) : 3 4 5 6 1 2 

SENSITIVITY  MEASURE  (1):  ALL  P ON  ONE  V=  0.560321229 

SENSITIVITY  MEASURE  (2):  ALL  P ON  ONE  V=  0.317454394 


VARIABLE  V=4  FUNCTION  VALUE  = 0.305936E+01 


TIME 

0.50 


IP 

1 


1 

LDVDP-1 
0. 213598E+00 
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2 

DVDP-2 

-0 . 703300E-02 


0.50 

2 

0 . 234289E+00 

-0.362882E-01 

0.50 

3 

-0 . 722285E-01 

0. 185719E+00 

0 . 50 

4 

-0. 112500E+00 

0 . 146042E+01 

0.50 

5 

-0. 110391E+00 

0.29182 1E-0 1 

0.50 

6 

-0. 164008E+00 

0 . 245461E+00 

PARAMETER  EFFECT  ORDER  (1): 
PARAMETER  EFFECT  ORDER  (2) : 


2 1 6 4 5 3 

6 4 3 2 5 1 


SENSITIVITY 

MEASURE 

(1)  : 

ALL 

P 

ON 

ONE 

V=  0.396828617 

SENSITIVITY 

MEASURE 

(2)  : 

ALL 

P 

ON 

ONE 

V=  1.493246967 

VARIABLE  V=5 

FUNCTION  VALUE  =-0 

. 428256E+01 

TIME  IP 

0.50  1 

0.50  2 

0.50  3 

0.50  4 

0.50  5 

0.50  6 

1 

LDVDP-1 
-0. 656031E+00 
0. 103553E+00 
0. 659698E+00 
-0.871288E-01 
0 . OOOOOOE+OO 
0. 000000E+00 

2 

DVDP-2 

0 . 688534E-01 
0. 15098  IE-01 
0.2582  3 3E+0 1 
-0. 176348E+01 
0. 000000E+00 
0. 000000E+00 

PARAMETER  EFFECT  ORDER  (1) : 3 

PARAMETER  EFFECT  ORDER  ( 2 ) : 3 

1 2 4 6 5 

4 12  6 5 

SENSITIVITY 

SENSITIVITY 

MEASURE  (1) : ALL  P 
MEASURE  (2) : ALL  P 

ON  ONE  V=  0.940155383 
ON  ONE  V=  2.370358344 

VARIABLE 

II 

> 

FUNCTION  VALUE  = 

0. 707297E+01 

1 

2 

TIME 

IP 

LDVDP-1 

DVDP-2 

0.50 

1 

0 . 220780E+00 

-0. 831865E-01 

0.50 

2 

0. 263169E+00 

-0.454338E-01 

0 .50 

3 

0 . 3 08584E+00 

-0. 238722E+01 

0 .50 

4 

-0.598843E-01 

0 . 204447E+01 

0.50 

5 

0. 100000E+01 

0. 313934E-12 

0 . 50 

6 

0. 000000E+00 

0. 000000E+00 

PARAMETER  EFFECT  ORDER  (1): 
PARAMETER  EFFECT  ORDER  (2): 


5 3 2 1 4 6 

3 4 1 2 5 6 


SENSITIVITY  MEASURE  (1): 
SENSITIVITY  MEASURE  (2): 


ALL  P ON  ONE  V=  1.103091942 
ALL  P ON  ONE  V=  3.526807997 
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2.2.5  Changing  Dynamic  Models 


Five  different  model  subroutines  are  included  on  disk  1 and  are 
stored  under  different  filenames.  Brief  descriptions  of  these 
models  are  given  below: 


SMD2 . FOR  - A 2 degree  of  freedom  spring-mass-damper  system 
used  to  demonstrate  model  development  for  primary  variables 
and  primary  parameters. 

SMD3 . FOR  - A 2 degree  of  freedom  spring-mass-damper  system 
used  to  demonstrate  model  development  for  secondary  variab- 
les and  primary  parameters. 

SMD4 . FOR  - A 2 degree  of  freedom  spring-mass-damper  system 
used  to  demonstrate  model  development  for  secondary  variab- 
les and  secondary  parameters. 

MODEL. FOR  - A 3 degree  of  freedom  nonlinear  vehicle  model 
based  on  vehicle  dynamics  simulations  developed  and 
copyrighted  by  Dr.  A.  G.  Nalecz^. 

ROLOVER. FOR  - An  eight  degree  of  freedom  vehicle  rollover 
model  based  on  the  simulation  developed  by  Systems 
Technology  Inc. 


The  particular  model  under  investigation  can  be  changed  by  compi- 
ling the  file  containing  the  model  desired  and  placing  it  in  the 
object  library  SENS.OLB  using  the  following  VAX  VMS  command: 


$ @LIB  TRNS  SMD3 


The  command  shown  above  will  compile  the  spring-mass-damper  model 
(File  SMD3.F0R)  shown  above,  and  replace  object  code  currently 
residing  in  the  object  library  SENS. 


Reprinted  from  "Lateral  Weight  Transfer  Simulation  (LWTS) " 
and  "Analysis  of  Suspension  Mechanics  Simulation  (ASMS) " by 
Dr.  A. G . Nalecz  by  permission  of  copyright  owner  (Dr.  A . G . 
Nalecz) . Year  of  first  publication  1986. 
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After  compilation  has  been  completed  the  user  must  once  again 
link  the  program  to  obtain  an  executable  version  of  the  sensiti- 
vity program.  This  is  accomplished  using  the  following  command. 


$ LINK  TRANS , TRNS/LIB, IMSLIBD/LIB 


Each  model  has  a parameter  data  file  associated  with  it.  The 
filenames  of  these  data  files  are  identical  to  the  filenames  of 
the  model  subroutines. 


Model  File  Name 


Data  File  Name 


MODEL. FOR 
RO LOVER. FOR 
SMD2 . FOR 
SMD3 . FOR 
SMD4 . FOR 


MODEL. DAT 
RO LOVER . DAT 
SMD2 . DAT 
SMD3 . DAT 
SMD4 . DAT 


These  data  files  include  the  parameter  values,  and  also  indicate 
which  parameters  and  variables  will  be  included  in  the  sensitiv- 
ity analysis.  During  program  execution  the  sensitivity  algorithm 
looks  for  this  file  under  the  name  VDAT.DAT.  Therefore  the  data 
file  associated  with  the  model  must  be  renamed  to  VDAT.DAT  before 
program  execution.  This  can  be  done  using  the  copy  command.  An 
example  is  shown  below: 


$ COPY  SMD3.DAT  VDAT.DAT 
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3.  SENSITIVITY  ANALYSIS  OF  STEADY  STATE  MODELS 


****  STEDY . FOR  **** 


3 . 1 Program  Description 


Program  STEDY. FOR  calculates  first  order  standard,  percentage  and 
logarithmic  sensitivity  functions  as  well  as  second  order 
sensitivity  functions  of  steady  state  vehicle  models.  The 
program  computes  the  parametric  sensitivity  functions  defined  as 
the  partial  derivatives  of  system  variables  with  respect  to 
system  parameters.  In  steady  state  analysis  the  system  must  take 
the  mathematical  form  of  a set  of  simultaneous  equations  (either 
linear  or  nonlinear) . The  algorithm  calculates  sensitivity 
functions  over  a range  of  wheel  steer  angles. 

Like  the  transient  sensitivity  algorithm  the  steady  state  sens- 
itivity algorithm  stores  its  results  in  a user  readable  printer 
output  file  as  well  as  a plotter  output  file.  First  order 
standard,  percentage,  logarithmic  and  second  order  standard 
sensitivity  functions  are  stored  in  each  file. 

The  main  program  STEDY. FOR  uses  several  subroutines  including  the 
IMSL  routine  ZSPOW.  The  same  BASIC  program  used  to  plot  output 
for  the  transient  sensitivity  algorithm  can  be  used  to  plot 
output  from  STEDY. FOR.  Short  descriptions  of  each  subroutine  are 
given  below. 


ZSPOW. OBJ  - An  IMSL  subroutine  used  to  solve  sets  of 
simultaneous  linear  and  nonlinear  equations.  The  routine  is 
based  on  Powell's  method. 

SSFCN . FOR  - A function  subroutine  required  by  IMSL  routine 
ZSPOW. 

INPUT2 . FOR  - Provides  interactive  input  for  data  required  by 
the  sensitivity  algorithm.  This  information  includes  data 
required  by  the  IMSL  routine  ZSPOW,  the  increment  of  steer 
angle  for  printer  and  plotter  output  storage.  The  output 
filenames,  the  range  of  steer  angle  desired,  and  the  user's 
initial  guess  for  a solution  to  the  set  of  simultaneous 
equations  are  also  specified  in  this  subroutine. 

VDATA . FOR  - Reads  the  vehicle  data  from  a user  generated 
data  file.  This  file  includes  the  values  for  all  primary 
vehicle  parameters,  and  also  indicates  which  variables  and 
parameters  (primary  and  secondary)  the  user  wishes  to 
include  in  sensitivity  calculations. 

SSMODEL.FOR  - Steady  state  three  degree  of  freedom  vehicle 
model.  This  subroutine  calculates  the  simultaneous 


31 


nonlinear  equations  which  are  solved  by  IMSL  routine  ZSPOW. 
The  subroutine  also  calculates  secondary  parameters  and 
secondary  variables  needed  to  compute  sensitivity  functions. 
SSMODEL  is  interchangeable,  allowing  the  sensitivity 
algorithm  to  analyze  various  systems. 


3 . 2 Steady  State  Model  Development 


In  general,  the  equations  of  motion  of  any  mechanical  system  can 
be  written  as  a system  of  first  ordesr  differential  equations  in 
the  following  form: 


— f -j_  ( x^,  X2  / » • « , x^ , t ) 


X2  f- 2 ( X2  i . « » , x^ , t ) 


(1) 


'n  fj-|  ( Xlf  X2  / • • • / XR,  t ) 


where  X(XpX2,...,xn)  represents  the  vector  of  state  variables 
which  describe  the  dynamic  state  of  the  model.  These  equations 
can  be  integrated  to  obtain  the  time  response  of  the  mechanical 
system.  In  steady  state  the  left  sides  of  these  equations  become 
equal  to  zero,  resulting  in  the  following  set  of  simultaneous 
algebraic  equations. 


0 = fx 

( Xx, 

X2  * • • • r Xj^  , t 

) 

0 = f2 

( xx. 

X2  / • • • / Xj^  , t 

) 

(2) 

0 = fn  ( xx,  x2, . • • , xn,  t ) 


These  algebraic  equations  can  be  either  linear  or  nonlinear 
depending  on  the  structure  of  the  mechanical  system  and  its 
mathematical  representation. 
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The  major  difference  between  analysis  of  a mechanical  system  in 
time  domain  and  analysis  of  steady  state  is  that  the  system 
transient  response  is  found  by  integrating  a system  of  simul- 
taneous differential  equations  and  the  system  steady  state 
response  is  obtained  by  solving  a set  of  simultaneous  algebraic 
equations.  Once  the  equations  of  motion  (1)  of  the  mechanical 
system  have  been  brought  to  form  (2)  it  is  possible  to  write  the 
steady  state  model  subroutine  required  by  the  sensitivity 
analysis . 

The  goal  of  a transient  vehicle  model  subroutine  is  to  determine 
the  time  rates  of  change  of  the  state  variables  (see  chapter  2.2) 
while  the  goal  of  the  steady  state  vehicle  model  subroutine  is  to 
calculate  the  right  hand  sides  of  the  equations  shown  above.  The 
technical  report  includes  listings  of  both  transient  and  steady 
state  model  subroutines  for  the  3 degree  of  freedom  vehicle  model 
analyzed  in  this  project.  These  listings  show  that  there  are 
very  few  differences  between  the  subroutine  structure  of  the 
transient  sensitivity  model  and  that  of  the  steady  state  model. 

The  variables  passed  to  subroutine  SSMODEL. FOR  through  the 
argument  list  are  explained  below: 


F - The  value  computed  for  the  right  hand  side  of  the 
steady  state  equations  of  motion  based  on  the  current 
estimate  of  state  vector  X. 


X - The  program's  solution  vector  of  the  simultaneous 
algebraic  equations. 

VAR  - The  vector  of  primary  and  secondary  variables. 

PAR  - A vector  of  model  inputs.  The  current  front  wheel 
steer  angle  resides  in  element  1 of  this  vector. 

FLG1  - An  dummy  integer  flag  which  is  used  for  debugging 
purposes . 


The  current  version  of  the  steady  state  sensitivity  program  is 
capable  of  analyzing  systems  which  are  represented  by  no  more 
than  10  simultaneous  non-linear  equations  of  motion.  The  program 
can  examine  the  sensitivity  of  up  to  50  variables.  Because  of 
these  limitations  the  model  must  use  the  following  array  sizes 
when  developing  the  model: 
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X(10) 

F ( 10) 
VAR (50) 
PAR (10) 


If  the  user  wishes  to  analyze  systems  which  exceed  the 
limitations  listed  above  then  he  must  redimension  arrays  X,  F, 
VAR,  and  PAR  in  the  following  subroutines  to  account  for  the 
increased  model  size. 


SSMODEL. FOR 
STEDY . FOR 
SSFCN . FOR 
INPUT2 . FOR 


It  should  be  noted  that  IMSL  routine  ZSPOW  does  not  guarantee  a 
solution  to  sets  of  nonlinear  simultaneous  equations  and  the 
chances  of  success  are  dependent  on  the  size  of  the  problem,  the 
nonlinearity  of  the  equations,  and  the  accuracy  of  the  initial 
guess  supplied  by  the  user. 


The  steady  state  sensitivity  program  defines  primary  variables  as 
state  variables  obtained  directly  from  the  solution  of  the 
simultaneous  algebraic  equations.  Secondary  variables  are 
defined  as  variables  which  are  functions  of  primary  variables  and 
primary  parameters.  In  the  steady  state  model  subroutine  the 
secondary  variables  are  declared  in  the  same  manner  as  used  in 
the  transient  model  subroutine. 

Primary  (constant)  parameters  and  secondary  (non-constant) 
parameters  are  also  treated  in  the  same  manner  used  in  the 
transient  model  subroutine.  Values  of  primary  parameters  and 
secondary  parameter  coefficients  are  passed  to  the  model  sub- 
routine through  the  common  block  VDAT.  The  values  of  secondary 
parameters  are  passed  back  to  the  sensitivity  algorithm  through 
common  block  SECP.  As  before  the  order  of  the  parameters  in 
common  block  VDATA  must  correspond  to  the  order  in  which  they 
appear  in  parameter  data  file  VDAT . DAT . Also  the  order  of 
secondary  parameters  in  common  block  SECP  must  correspond  to  the 
order  of  appearance  of  secondary  parameter  coefficients  in  common 
block  VDATA.  After  each  secondary  parameter  has  been  computed  in 
the  model  subroutine  it  should  be  multiplied  by  the  secondary 
parameter  coefficient  associated  with  it  (see  chapter  2.2.4). 

The  structure  of  the  parameter  data  file  VDAT . DAT  used  in  the 
steady  state  sensitivity  analysis  is  quite  similar  to  that  used 
in  the  sensitivity  analysis  in  the  time  domain.  This  structure 
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is  described  in  the  instructions  for  the  transient  sensitivity 
algorithms  (chapter  2.2).  An  example  of  the  vehicle  data  file 
VDAT.DAT  used  for  the  three  degree  of  freedom  vehicle  model  is 
shown  below: 


1/ 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

1, 

0, 

0, 

1, 

0, 

If 

1 

/ 

If 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


DIST  FROM  CG  TO  FRONT  AXLE 

FORWARD  ACCELERATION  (+  BRAKING) 

DIST  FROM  CG  TO  REAR  AXLE 

TOTAL  SUSPENSION  DAMPING 

FRONT  WHEEL  STEERING  ANGLE 

REAR  WHEEL  STEERING  ANGLE 

STEER  OF  REAR  WHEEL/BODY  ROLL 

FRONT  ROLL  CENTER  HEIGHT 

RATE  OF  FRONT  WHEEL  INCL.  W/  ROLL 

CG  HEIGHT  OF  SPRUNG  MASS 

CG  HEIGHT  OF  FRONT  UNSPNG  MASS 

CG  HEIGHT  OF  REAR  UNSPNG  MASS 

DIST  FROM  ROLL  AXIS  TO  CG 

ROLL  MOMENT  OF  INERTIA  IX 

YAW  MOMENT  OF  INERTIA  IZ 

FRONT  SUSPENSION  ROLL  STIFFNESS 

REAR  SUSPENSION  ROLL  STIFFNESS 

TOTAL  SUSPENSION  ROLL  STIFFNESS 

WHEEL  BASE 

VEHICLE  MASS 

VEHICLE  SPRUNG  MASS 

VEHICLE  FRONT  UNSPRUNG  MASS 

VEHICLE  REAR  UNSPRUNG  MASS 

FRONT/REAR  BRAKING  PROPORTION 

HEAVY  BRAKING  PROPORTION 

REAR  ROLL  CENTER  HEIGHT 

FRONT  TRACK  WIDTH 

REAR  TRACK  WIDTH 

FORWARD  SPEED 

FRONT  WHEEL  DRIVE  (1)  REAR  (2) 
CALSPAN  COEFFICIENT  A0 
CALSPAN  COEFFICIENT  A1 
CALSPAN  COEFFICIENT  A2 
CALSPAN  COEFFICIENT  B1 
CALSPAN  COEFFICIENT  B3 
CALSPAN  COEFFICIENT  B4 
2 X FRONT  ALIGNING  STIFFNESS 
2 X REAR  ALIGNING  STIFFNESS 
INCLINATION  STIFFNESS 
CALSPAN  COEFFICIENT  P0 
CALSPAN  COEFFICIENT  PI 
CALSPAN  COEFFICIENT  P2 
SKID  NUMBER 

CONVERSION  FACTOR  FOR  FORCE 
GRAVITATIONAL  ACCELERATION 
GRAVITATIONAL  CONSTANT 
ENDPP 


' , 1.2D0 
' , 0 . 0D0 
' , 1. 3D0 
' , 37500. 0D0 
' , 0 . 01745D0 
' , 0. 0D0 
' , 0 . 05D0 
' , 0. 1711D0 
' , 0. 1D0 
' , 0. 6D0 
' , 0. 3D0 
' , 0. 3D0 
' , 0.5099D0 
' , 1100 . 0D0 
' , 1300. 0D0 
' , 50000. 0D0 
' , 30000. 0D0 
' , 80000. 0D0 
' , 2 . 5D0 
' , 1000. 0D0 
' , 900 . 0D0 
' , 60. 0D0 
' ,40. 0D0 
' ,0.6D0 
' , 0. 1D0 
' , 0. 01D0 
' , 1. 3D0 
' , 1. 3D0 
' , 25. 0D0 
' , 2 . 0D0 
' , 2380. 2D0 
' , 9 . 36D0 
' , 3621. 2D0 
' , -2 . 876D-04 
' , 1.274D0 
' , 2 . 6550D-08 
' , 3000. 0D0 
' , 3000. 0D0 
' , 6000. 0D0 
' , 1. 1972D0 
• , -2 . 992D-04 
' , 5.205D-08 
' , 80. 0D0 
' , 4 . 448D0 
' , 9 . 81D0 
' , 1. 0D0 
' , 0. 0D0 
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0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


PCF1  CORNERING  STIFFNESS 

1 

1 

, 1. 0D0 

PCF2  CORNERING  STIFFNESS 

2 

8 

, 1. 0D0 

PCF3  CORNERING  STIFFNESS 

3 

I 

, 1. 0D0 

PCF4  CORNERING  STIFFNESS 

4 

8 

, 1. ODO 

PFY1  FY1 

8 

, 1. ODO 

PFY2  FY2 

8 

, 1. 0D0 

PFY3  FY3 

8 

, 1. ODO 

PFY4  FY4 

0 

, 1. ODO 

PAB1  ALPHA  BAR  1 

8 

, 1. ODO 

PAB2  ALPHA  BAR  2 

8 

, 1. ODO 

PAB3  ALPHA  BAR  3 

8 

, 1. ODO 

PAB4  ALPHA  BAR  4 

8 

, 1. ODO 

PFA1  F OF  ALPHA  BAR 

1 

8 

, 1 . ODO 

PFA2  F OF  ALPHA  BAR 

2 

8 

, 1. ODO 

PFA3  F OF  ALPHA  BAR 

3 

8 

, 1. ODO 

PFA4  F OF  ALPHA  BAR 

4 

8 

, 1. ODO 

PMUX1  X FRICTION  COEFF. 

1 

8 

, 1. ODO 

PMUX2  X FRICTION  COEFF. 

2 

8 

, 1. ODO 

PMUX3  X FRICTION  COEFF. 

3 

8 

, 1. ODO 

PMUX4  X FRICTION  COEFF. 

4 

8 

, 1. ODO 

PMUY1  Y FRICTION  COEFF. 

1 

8 

, 1. ODO 

PMUY2  Y FRICTION  COEFF. 

2 

8 

, 1. ODO 

PMUY3  Y FRICTION  COEFF. 

3 

8 

, 1. ODO 

PMUY4  Y FRICTION  COEFF. 

4 

8 

, 1. ODO 

PWTF  FRONT  LATERAL 

WT. 

TRANSFER 

8 

, 1. ODO 

PWTR  REAR  LATERAL  WT.  TRANSFER 

f 

, 1. ODO 

ENDSP 

8 

, 0. ODO 

1 V 

8 

, 0. ODO 

2 r 

0 

, 0. ODO 

3 FAI 

0 

, 0. ODO 

4 FRONT  SLIP  ANGLE 

8 

, 0 . ODO 

5 REAR  SLIP  ANGLE 

8 

, 0. ODO 

6 FY1 

8 

, 0. ODO 

7 FY2 

0 

, 0 . ODO 

8 FY3 

0 

, 0 . ODO 

9 FY4 

8 

, 0. ODO 

10  SMAX1 

0 

, 0. ODO 

11  SMAX2 

0 

, 0. ODO 

12  SMAX3 

8 

, 0. ODO 

13  SMAX4 

0 

, 0. ODO 

ENDOV 

0 

, 0. ODO 

3 . 3 Execution  of  Steady  State  Sensitivity  Program 


When  the  user  has  written  his  steady  state  model  subroutine 
SSMODEL. FOR  and  generated  the  parameter  data  file  VDAT.DAT  he 
must  compile  the  model  and  link  the  subroutines  together  to  form 
an  executable  file  of  the  steady  state  sensitivity  algorithm. 
This  can  be  done  using  the  following  commands: 
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$ @LIB  SDST  SSMODEL 


$ LINK  STEDY, SDST/LIB, IMSLIBD/LIB 


Once  the  executable  version  is  obtained  the  program  can  be  run 
using  the  following  command: 


$ RUN  STEDY 


***  Data  Entry  *** 

1 . Enter  From  Keyboard 

2.  Read  From  File 

1 (For  input  from  keyboard) 


As  before  the  user  has  the  choice  of  entering  simulation  data 
from  the  keyboard  or  reading  it  from  an  existing  data  file.  For 
this  example  the  data  will  be  entered  using  the  keyboard. 


Input  Number  of  Equations  to  Solve 
2 


This  represents  the  number  of  simultaneous  algebraic  equations 
which  must  be  solved. 


Initial  Solution  Estimate 

X(  1)  = 0.0D0 
X(  2)  = 0.0D0 


Do  you  wish  to  change  these  parameters  (Y/N) 
N 


The  user  must  input  an  initial  guess  for  the  solution  to  the 
simultaneous  algebraic  equations.  The  default  values  of  zero 
will  be  used  for  this  example. 
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***  Data  Required  By  IMSL  Routine  ZSPOW  **** 

Default  Values 

NSIG  =4 
ITMAX=2  00 

Do  you  wish  to  change  any  of  these  parameters  (Y/N) 
N 


These  parameters  are  required  by  IMSL  routine  ZSPOW.  NSIG 
specifies  the  accuracy  level  of  the  routine.  ITMAX  represents  the 
maximum  allowable  number  of  function  evaluations  which  can  be 
performed  during  each  call  to  ZSPOW.  The  user  of  the  sensitivity 
algorithm  has  the  option  of  accepting  these  default  values  or 
changing  them.  They  will  not  be  changed  for  this  example.  For 
more  information  regarding  these  parameters  please  consult  your 
IMSL  manual. 


Input  the  Number  of  State  Variables 
13 


This  represents  the  total  number  of  variables  (both  primary  and 
secondary)  which  are  to  be  included  in  the  sensitivity 
computations . 


***  Program  Parameters  *** 

Derivative  Method  - Central  Difference  With  3 Points 
Percent  used  in  Central  Difference  - 1.0D-2 


Do  You  Wish  to  Change  These  Parameters  (Y/N) 
N 


Here  the  user  has  the  option  of  accepting  the  default  values 
specifying  the  number  of  points  used  for  the  derivative 
evaluation  and  the  percentage  parameter  change  associated  with 
sensitivity  function  calculation.  The  user  may  change  these 
values  if  he  so  desires.  The  values  will  not  be  changed  for  this 
example. 
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***  Final  Steer  Angle  and  Steer  Increment  *** 


Input  Total  Steering  Angle  in  Degrees  (Double  Precision) 

1.5D0 

Input  Output  Steering  Increment  in  Degrees  (Double  Precision) 

0. 3D0 

Input  Filename  For  Printer  Output  Data 
STEDY . OUT 

Would  You  Like  the  Program  to  Create  a Plotting  Output  File  (Y/N) 
Y 

Input  Storage  time  Increment  for  Plot  (Double  Precision) 

0.1D0 

Input  Filename  For  Plotting  Data 
STEDY . PLT 


The  inputs  shown  above  indicate  that  the  sensitivity  functions 
will  be  calculated  for  a front  wheel  steer  angle  range  of  0 to 
1.5  degrees.  The  sensitivity  functions  will  be  written  to  the 
printer  output  data  file  every  0.3  degrees  of  front  wheel 
steering  input  and  the  printer  output  data  filename  will  be 
STEDY. OUT.  A plotter  file  is  to  be  created  and  plotting  data 
will  be  stored  under  the  filename  STEDY. PLT.  The  sensitivity 
functions  will  be  stored  in  the  plotting  output  file  at  every  0.1 
degrees  of  front  wheel  steering  input. 


Do  You  Want  to  Store  This  Data  in  a File  (Y/N) 

Y 

What  Filename  Do  You  Wish  to  Give  the  Input  Data 
STEDY . INI 


This  allows  the  user  to  store  the  information  entered  above  into 
a data  file.  Once  the  data  is  stored  in  a file  the  user  can  let 
the  program  read  the  data  from  the  file  on  subsequent  runs, 
rather  than  inputing  the  data  from  the  keyboard  each  time  (see 
data  entry  choice) . In  this  case  the  initial  data  file  was 
stored  in  a file  called  STEDY . INI . 

Sample  results  obtained  from  the  three  degree  of  freedom  steady 
state  vehicle  model  (written  in  file  STEDY. OUT)  after  the  first 
printing  steer  increment  are  shown  below: 
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STATE  VARIABLE  V 


1,  FUNCTION  VALUE  =-0 . 290125E+00 


1 

2 

STEER 

IP 

DV-1 

DVDP-1 

0.300 

1 

-0. 115850E-02 

-0 . 965418E-01 

0.300 

2 

-0 . 961352E-08 

”0. 160225E-05 

0.300 

3 

-0.348288E-03 

-0. 683053E-01 

0.300 

4 

“0 . 9613  52E-08 

-0 . 7 3 950 IE- 09 

0.300 

5 

0 . 569249E-04 

0. 113850E-06 

0.300 

6 

”0. 285778E-04 

-0. 952593E-07 

3 

LDVDP-1 
0. 399311E+00 
0 . 331358E-05 
0. 120048E+00 
0. 331358E-05 
-0. 196208E-01 
0. 985017E-02 


4 

DVDP-2 

0 . 301453E+00 
0. 534084E-03 
-0.249900E-01 
0. 113769E-09 
0. 154541E-11 
-0.284568E-11 


PARAMETER  EFFECT  ORDER  (1):  1 3 
PARAMETER  EFFECT  ORDER  (2):  1 3 
PARAMETER  EFFECT  ORDER  ( 3 ) : 1 3 
PARAMETER  EFFECT  ORDER  (4) : 1 3 


5 6 4 2 
2 5 6 4 
5 6 2 4 
2 4 6 5 


SENSITIVITY  MEASURE 
SENSITIVITY  MEASURE 
SENSITIVITY  MEASURE 
SENSITIVITY  MEASURE 


(1)  : 

ALL 

P 

ON 

ONE 

(2)  : 

ALL 

P 

ON 

ONE 

(3)  : 

ALL 

P 

ON 

ONE 

(4)  s 

ALL 

P 

ON 

ONE 

V=  1 . 20924297E-03 
V=  1. 18260665E-01 
V=  4 o 17427906E-01 
V=  7. 13390714E-01 


The  format  used  to  present  these  results  is  identical  to  that 
used  for  the  transient  sensitivity  analysis  (see  chapter  2.2.1). 
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4.  SENSITIVITY  ANALYSIS  IN  THE  FREQUENCY  DOMAIN 


****  FREQDOM3 . FOR  **** 


4 . 1 Program  Description 


Program  FREQDOM3 . FOR  calculates  first  order  standard,  percentage 
and  logarithmic  sensitivity  functions  as  well  as  second  order 
standard  and  logarithmic  sensitivity  functions  of  linear  vehicle 
models  written  in  the  frequency  domain  over  a user  specified 
frequency  range.  The  program  computes  the  parametric  sensitivity 
functions  defined  as  the  partial  derivative  of  elements  of  the 
system's  transfer  function  matrix  with  respect  to  system 
parameters.  The  system  must  be  written  as  a linear  set  of 
differential  equations. 

The  program  stores  its  results  in  two  data  files.  The  first  file 
is  a user  readable  printer  output  data  file.  The  values  of  user 
selected  elements  of  the  transfer  function  matrix  are  stored  in 
this  file  along  with  standard,  percentage,  logarithmic 

first  order  sensitivity  functions  as  well  as  standard  and 
logarithmic  second  order  sensitivity  functions.  The  file  also 
shows  general  sensitivity  measures  and  the  order  of  influence  of 
system  parameters  on  elements  of  the  transfer  function  matrix 
based  on  each  type  of  sensitivity  function.  Hardcopy  of  this 
file  can  be  produced  by  sending  it  to  a line  printer.  The  second 
data  file  is  a plotter  output  data  file.  This  file  contains 
values  for  selected  elements  of  the  transfer  function  matrix  and 
values  of  first  order  standard,  percentage  and  logarithmic 
sensitivity  functions  as  well  as  second  order  standard 
sensitivity  functions.  Output  data  is  usually  stored  more  fre- 
quently in  the  plotter  output  data  file  than  in  the  printer  data 
output  file  (using  a smaller  plotting  increment  step  for  better 
graph  resolution)  and  is  stored  using  a file  structure  which  is 
more  compact  than  that  used  for  the  printer  output  data  file. 
The  same  BASIC  program  used  to  plot  output  in  the  previous  sensi- 
tivity algorithms  can  be  used  to  plot  output  from  FREQD0M3 . FOR. 

The  frequency  domain  sensitivity  algorithm  is  restricted  to 
analysis  of  primary  parameters  (parameters  which  remain  constant) 
and  primary  variables  (variables  which  are  elements  of  the 
transfer  function  matrix) . 

The  main  program  FREQD0M3 . FOR  uses  several  subroutines;  short 
descriptions  of  each  subroutine  are  given  below. 


INPUT 3 . FOR  - Provides  interactive  input  for  the  data  re- 
quired by  the  sensitivity  algorithm.  This  information 
includes  the  frequency  range  over  which  the  sensitivity 
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functions  will  be  computed,  the  frequency  increments  of 
printer  and  plotter  data  storage,  the  data  output  file 
names . 

FVDATA . FOR  - Reads  the  vehicle  parameters  from  a user 
generated  data  file.  This  file  includes  the  values  for  all 
primary  vehicle  parameters,  and  also  flags  which  elements  of 
the  transfer  function  matrix  and  parameters  the  user  wishes 
to  include  in  the  sensitivity  computations. 

FMODEL. FOR  - Calculates  the  coefficient  matrix  [A]  and  the 
input  matrix  [B]  of  the  linear  vehicle  model  written  in  the 
form  [X]  = [A] [X]  + [B] [C]  where  [A]  is  a square  matrix,  [X] 
is  a vector  state  variables,  [B]  is  a column  matrix,  and  [C] 
represents  the  system  inputs. 

TRANSFER. FOR  - Computes  the  complex  transfer  function  matrix 
given  the  coefficient  matrix  [A] , the  input  matrix  [B] , and 
the  frequency  of  input  excitation. 


4 . 2 Frequency  Domain  Sensitivity  Model  Development 


In  order  to  analyze  the  sensitivity  of  a model  written  in  the 
frequency  domain  the  user  must  develop  a model  according  to 
certain  guidelines  and  must  also  create  a parameter  data  file 
using  a set  structure. 

In  order  to  analyze  the  sensitivity  of  a mechanical  system  in  the 
frequency  domain  the  system  must  be  written  as  a set  of  linear 
first  order  differential  equations: 


[X]  = [A] [ X]  + [B][C] 


Where  [X]  represents  the  vector  of  state  variables,  [A]  and  [B] 
represent  matricies  whose  elements  are  functions  of  the  system 
parameters,  and  [C]  is  a vector  which  represents  the  input  given 
to  the  system  (e.g.  steering  angle) . Once  we  are  able  to  calc- 
ulate matrices  [A]  and  [B]  we  are  able  to  calculate  the  system's 
transfer  function  matrix.  The  transfer  function  matrix  is  used 
to  obtain  the  amplification  factor  and  phase  angle  of  the  system 
based  on  the  frequency  of  input  excitation.  The  sensitivity 
algorithm  computes  the  change  of  elements  of  the  complex  transfer 
function  matrix  caused  by  small  parameter  changes. 

The  subroutine  listing  of  a three  degree  of  freedom  vehicle  model 
is  shown  below: 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


★ •A***************'******'*****'**'*'*'*****'***'**'*'*'**'*'*'*'*'***'**** 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


SUBROUTINE  - FMODEL 


* 
* 
* 

PURPOSE  - To  return  coefficient  matricies  [A]  & [B]  * 

of  3 DOF  linear  steady  state  vehicle  model  * 
of  the  form:  * 

. * 

[X]=[A] [X]+[B] [C]  * 

* 

REQUIRES  DATA  - FROM  COMMON  BLOCK  VDAT  * 


C 


c 


c 

c 

c 

c 


c 

c 


c 


c 


SUBROUTINE  FMODEL (N, A, B) 

REAL  A ( 10, 10) , B (10, 10) , L, M , MS , K13 , K2 3 , K3 3 , KU1 , KU2 , IZ , 
IX , K, Mil , Ml 3 , M2 2 , M3 1 , M3  3 , Mill , MI 13 , MI 2 2 , MI 3 1 , MI 3 3 


COMMON/VDAT/ 

Al, 

CTOT , 

Dl, 

DF, 

DI, 

H, 

• K, 

IX, 

IZ, 

L, 

M, 

MS, 

u. 

CMF, 

! CMR, 

CT, 

CA, 

CB 

CALL  FVDATA 

G=9 . 8 1 

B1=L-A1 

M11=M 
M13=MS*H 
M22=IZ 
M3 1=MS*H 
M3  3=IX 

DET=M11*M22*M33-M13*M22*M31 

MI 1 1=M2  2 *M3  3 / DET 

MI13=-M13*M22/DET 

MI 2 2= (M11*M33-M13*M31) /DET 

MI31=-M22*M31/DET 

MI 3 3 =M1 1 *M2  2 / DET 

C11=CA*M*G* (Al+Bl) /U/L+4 . 0*CB/U 
C12=2 . 0*CB* (Al-Bl) /U+M*U 

C2 1=2 . 0*CB*Al/U-2 . 0*CB*B1/U+ ( CMF+CMR) /U 

C22=CA*Al*Al*M*G*Bl/U/L+2 . 0*A1*A1*CB/U+CA*B1*B1*M*G*A1/U  L- 
2 . 0*B1*B1*CA/U+A1*CMF/U-B1*CMR/U 
C32=MS*H*U 
C33=CTOT 
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o o o non 


K13=DI*CT+DF*M*G*A1*CA/L+DF*2 . 0*CB 

K23=DF*M*G*A1*B1*CA/L+DF*2 . 0*CB*B1-DI*CF*A1-DI*CMR 
K33=K-MS*G*H 

KUl=M*G*Bl*CA/L+2 . 0*CB 
KU2=M*G*B1*CA*A1+2*CB*A1+CMF 


GENERATE  MATRIX  [A] 


A ( 1 , 1) =-Cll*MIll 

A(1,2)=-C32*MI13 

A ( 1 , 3 ) =-C33*MI13 

A(1,4)=-K13*MI11-K33*MI13 

A(2, 1)=-C21*MI22 

A(2,2)=-C22*MI22 

A ( 2 , 3 ) =0 . 0 

A ( 2 , 4 ) =-K2  3 *MI22 

A ( 3 , 1)=-C11*MI31 

A ( 3 , 2 ) =-C12*MI31“C32*MI33 

A (3 , 3 ) =-C33*MI33 

A(3 , 4) =-K13*MI31-K33*MI33 

A(4, 1)=0.0 

A(4/2)=0c0 

A ( 4 , 3 ) =1 . 0 

A (4, 4) =0.0 


GENERATE  MATRIX  [B] 

B(l, 1)=KU1*MI11 
B (2 , 1) =KU2*MI22 
B ( 3 / 1) =KU1*MI3 1 
B ( 4 , 1 ) =0 . 0 


RETURN 

END 


The  variables  in  the  subroutine  argument  list  have  the  following 
meanings : 


N - The  number  of  first  order  linear  differential  equations 
describing  the  model. 

A - The  N x N matrix  [A]  in  the  expression  describing  the 
system's  differential  equations  [X]  = [A] [X]  + [B] [C] 

B - The  N x M vector  [B]  (where  M represents  the  number  of 
inputs  given  to  the  system)  in  the  expression  describing 
the  system's  differential  equations  [X]  = [A] [X]  + 

[B] [C] 
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As  before,  the  common  block  VDAT  is  used  to  pass  the  system 
parameters  into  the  subroutine.  The  order  in  which  the 
parameters  appear  in  common  block  VDAT.DAT  must  be  identical  to 
their  order  of  appearance  in  data  file  FVDAT.DAT. 

Matricies  [A]  and  [B]  must  be  dimensioned  10  x 10.  If  the  user 
wishes  to  develop  a model  which  is  described  by  more  than  10 
first  order  differential  equations,  or  has  more  than  10  inputs 
then  the  user  must  redimension  matricies  [A]  and  [B]  in  the  main 
program  FREQDOM3 . FOR  and  in  subroutine  TRANSFER. FOR  to 
accommodate  larger  models. 

The  first  executable  line  in  the  model  subroutine  should  be  a 
call  to  subroutine  FVDATA.FOR.  This  subroutine  reads  the  system 
parameters  from  data  file  FVDAT.DAT  and  passes  them  into  the 
model  subroutine  through  common  block  VDAT.DAT.  After  the  call 
to  subroutine  FVDATA.FOR  has  been  made  remainder  of  the  model 
subroutine  is  used  to  calculate  matricies  [A]  and  [B] . 

The  structure  of  vehicle  data  file  FVDAT.DAT  is  quite  similar  to 
that  used  in  the  transient  and  steady  state  sensitivity 
algorithms  (VDAT.DAT),  however  several  differences  exist.  The 
file  FVDAT.DAT  begins  with  the  number  of  first  order  differential 
equations  used  in  the  model;  for  the  example  model  this  number  is 
4 . The  next  number  represents  the  number  of  inputs  given  to  the 
model;  for  the  example  model  this  number  is  1.  After  these 
numbers  are  placed  at  the  top  of  the  vehicle  data  file  the 
primary  parameters  are  listed  using  a format  identical  to  that 
used  in  the  transient  and  steady  state  sensitivity  algorithms. 
The  parameter  description  ENDPP  is  again  used  to  note  the  end  of 
primary  parameters.  The  frequency  domain  sensitivity  algorithm 
is  not  capable  of  calculating  sensitivity  functions  with  respect 
to  secondary  parameters,  therefore  there  is  no  secondary 
parameter  flag  ENDSP  in  file  FVDAT.DAT. 

The  frequency  domain  sensitivity  algorithm  computes  influence  of 
system  parameters  on  one  or  more  elements  of  the  transfer 
function  matrix.  The  user  specifies  which  elements  of  the 
transfer  function  matrix  he  wishes  to  include  in  the  sensitivity 
calculations  by  listing  the  row  and  column  indices  of  each 
element  at  the  bottom  of  file  FVDAT.DAT.  The  final  entry  in  data 
file  FVDAT.DAT  is  row  and  column  indices  0,0.  This  entry  signals 
the  program  that  there  are  no  more  elements.  The  sample  data 
file  below  shows  that  the  sensitivity  functions  will  be  computed 
for  transfer  function  elements  (1,1)  and  (2,1)  with  respect  to 
those  parameters  having  a 1 in  their  first  data  field.  In 
general  transfer  function  element  i,j  represents  the  input/output 
relationship  of  state  variable  i to  system  input  j . For  example 
transfer  element  (2,1)  represents  the  input/ output  relationship 
between  the  state  variable  x2  and  system  input  1. 
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4 

1 

1, 

0, 

0, 

1/ 

0, 

1, 

0, 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


DIST  FROM  CG  TO  FRONT  AXLE 
TOTAL  SUSPENSION  ROLL  DAMPENING 
FRONT  WHEEL  STEERING  ANGLE 
STEER  OF  REAR  WHEEL/ BODY  ROLL 
RATE  OF  FRONT  WHEEL  INCL.  W/  ROLL 
HEIGHT  OF  SPRUNG  MASS 
TOTAL  SUSPENSION  ROLL  STIFFNESS 
IX 
IZ 

WHEEL  BASE 
VEHICLE  MASS 
VEHICLE  SPRUNG  MASS 
FORWARD  SPEED 

2 X FRONT  ALIGNING  STIFFNESS 
2 X REAR  ALIGNING  STIFFNESS 
INCLINATION  STIFFNESS 
CORNERING  STIFFNESS  COEFF.  A 
CORNERING  STIFFNESS  COEFF.  B 
ENDPP 


1,1 
2 , 1 
0,0 


,0.869 
, 1600.0 
,0.01745 
,0.02 
,0.2 
,0.5967 
, 34383 . 0 
,400.0 
,2000.90 
, 2 . 66 
, 1230.16 
, 1095.23 
,25.0 
, 3000 . 0 
,3000.0 
, 6000.0 
,1.675 
, 17295.7 
,0.0 


4 . 3 Execution  of  Frequency  Domain  Sensitivity  Programs 


Once  the  user  has  written  his  model  subroutine  and  generated  the 
parameter  data  file  FVDAT.DAT  he  must  compile  the  model  and  link 
the  subroutines  together  to  form  an  executable  file  of  the 
sensitivity  algorithm.  This  can  be  done  using  the  following 
commands : 


$ @LIB  FREQ  FMODEL 
$ LINK  FREQDOM3 , FREQ/ LIB 


Once  the  executable  version  is  obtained  the  program  can  be  run 
using  the  following  command: 


$ RUN  FREQDOM3 
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***  Data  Entry  *** 

1.  Enter  From  Keyboard 

2 . Read  From  File 

1 (For  input  from  keyboard) 

***  Sensitivity  Type  *** 

1.  Amplitude  Sensitivity 

2.  Phase  Sensitivity 

Input  the  Number  of  Your  Choice 
1 


If  the  user  chooses  to  calculate  amplitude  sensitivity  then  the 
program  will  calculate  the  sensitivity  functions  based  on  the 
magnitude  of  the  complex  element  of  the  transfer  function  matrix. 
If  the  user  selects  phase  sensitivity  then  the  program  will 
calculate  the  sensitivity  functions  based  on  the  phase  angle  of 
the  complex  element  of  the  transfer  function  matrix. 


***  Program  Parameters  *** 


Derivative  Method  - Central  Difference  With  3 Points 
Percent  used  in  Central  Difference  - 1.0E-2 

Do  You  Wish  to  Change  These  Parameters  (Y/N) 

N 

***  Frequency  Start,  Stop,  and  Step  *** 

Starting  Frequency  (Rad/S) 

0.0 

Ending  Frequency  (Rad/S) 

10.0 

Frequency  Increment  For  Output  Data  (Rad/S) 

2.50 

Input  Filename  For  Output  Data 
FREQ . OUT 

Would  You  Like  the  Program  to  Create  a Plotting  Output  File  (Y  N’) 
Y 

Input  Storage  Frequency  Increment  for  Plot  (Rad/S) 

0.50 
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Input  Filename  For  Plotting  Data 
FREQ . PLT 

Do  You  Want  to  Store  This  Data  in  a File  (Y/N) 

Y 

What  Filename  Do  You  Wish  to  Give  the  Input  Data 
FREQ . INI 


The  user  inputs  shown  above  indicate  that  the  sensitivity  algor- 
ithm will  calculate  sensitivity  functions  from  0 to  10  Rad/S  and 
will  store  printer  results  in  file  FREQ. OUT  every  2.5  Rad/S  and 
plotter  results  in  file  FREQ. PLT  every  0.5  Rad/S. 

Sample  output  from  file  FREQ. OUT  is  shown  below.  The  values  of 
the  second  order  logarithmic  sensitivity  function  are  omitted  due 
to  space  considerations. 


AMPLITUDE  OF  ELEMENT  C(  1,  1)  = 0.644184E+02 


1 

2 

3 

4 

FREQ 

IP 

DV-1 

DVDP-1 

LDVDP-1 

DVDP-2 

2 . 50 

1 

0. 197104E+01 

0.226817E+03 

0. 305974E+01 

0. 118064E+04 

2.50 

2 

0. 800209E-01 

0 . 400105E+03 

0. 124221E+00 

0.476837E+04 

2.50 

3 

0. 494783E+00 

0 . 829200E+02 

0 . 7 68077E+00 

-0. 641228E+03 

2 . 50 

4 

-0. 219418E+00 

-0. 548544E-01 

-0. 340613E+00 

-0. 191212E-03 

PARAMETER  EFFECT  ORDER 

(1)  : 

1 

3 

4 

2 

PARAMETER  EFFECT  ORDER 

(2)  : 

2 

1 

3 

4 

PARAMETER  EFFECT  ORDER 

(3)  : 

1 

3 

4 

2 

PARAMETER  EFFECT  ORDER 

(4)  : 

2 

1 

3 

4 

SENSITIVITY  MEASURE 

(1) 

: ALL 

P 

ON 

ONE 

V= 

2 . 046489 

SENSITIVITY  MEASURE 

(2) 

: ALL 

P 

ON 

ONE 

V= 

467.3228 

SENSITIVITY  MEASURE 

(3) 

: ALL 

P 

ON 

ONE 

v= 

3 . 166044 

SENSITIVITY  MEASURE 

(4) 

: ALL 

P 

ON 

ONE 

v= 

4956.817 

AMPLITUDE  OF  ELEMENT  C(  2,  1)  = 0.170322E+01 


1 

2 

3 

4 

FREQ 

IP 

DV-1 

DVDP-1 

LDVDP-1 

DVDP-2 

2 . 50 

1 

0 . 446776E-01 

0. 514126E+01 

0 . 262313E+01 

0. 268819E+02 

2 . 50 

2 

0. 206542E-02 

0. 103271E+02 

0. 121266E+00 

0 . 834465E+02 

2 . 50 

3 

-0. 244207E+00 

-0 . 4 092  63E+02 

-0. 143380E+02 

0. 807092E+02 

2 . 50 

4 

0. 122147E+00 

0. 305367E-01 

0. 717154E+01 

-0 . 696853E-04 
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PARAMETER  EFFECT  ORDER 

(1)  : 

3 

4 

1 

2 

PARAMETER  EFFECT  ORDER 

(2)  : 

3 

2 

1 

4 

PARAMETER  EFFECT  ORDER 

(3)  : 

3 

4 

1 

2 

PARAMETER  EFFECT  ORDER 

(4)  : 

2 

3 

1 

4 

SENSITIVITY  MEASURE 

(1) 

: ALL 

P 

ON 

ONE 

V= 

0.2766761 

SENSITIVITY  MEASURE 

(2) 

: ALL 

P 

ON 

ONE 

V= 

42.26542 

SENSITIVITY  MEASURE 

(3) 

: ALL 

P 

ON 

ONE 

V= 

16.24512 

SENSITIVITY  MEASURE 

(4) 

: ALL 

P 

ON 

ONE 

v= 

119 . 2226 

The  AMPLITUDE  OF  ELEMENT  C(l,l)  = 0.644184E+02  indicates  that  the 
sensitivity  results  which  follow  correspond  to  the  amplitude  of 
the  complex  transfer  function  element  0(1,1). 

The  numbers  listed  under  the  column  labeled  IP  indicate  which 
parameters  correspond  to  the  sensitivity  values.  IP=1  represents 
the  first  parameter  selected  in  data  file  VDAT.DAT  (with  the 
first  data  field  set  equal  to  1)  , IP=2  the  second,  IP=3  the 
third,  ecte 

The  columns  labeled  DV-l,  DVDP-1,  LDVDP-1,  and  DVDP-2  represent 
the  various  types  of  sensitivity  functions  of  variable  V with 
respect  to  parameter  IP.  A description  of  the  column  labels  are 
shown  below: 


DV-l 

DVDP-1 

LDVDP-1 

DVDP-2 

LDVDP-2 


First  Order  Percentage  Sensitivity 

First  Order  Standard  Sensitivity 

First  Order  Logarithmic  Sensitivity 

Second  Order  Standard  Sensitivity 

Second  Order  Logarithmic  Sensitivity  (Not  Shown) 


The  PARAMETER  EFFECT  ORDER  shows  the  order  of  influence  (from 
greatest  to  least)  of  all  parameters  on  variable  V based  on  each 
sensitivity  type.  The  SENSITIVITY  MEASURE  shows  the  combined 
effect  of  all  parameters  and  is  based  on  each  sensitivity  type. 
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5.  PLOTTING  SENSITIVITY  RESULTS 


5 . 1 Downloading  Plotting  Output  Data  File  To  IBM  PC  Using  Procomm 


(Note:  It  is  not  mandatory  that  operators  use  Procomm  communica- 
tion software  to  download  the  plotting  data  files.  Users  may  use 
another  communications  software  package  if  they  so  desire.) 

In  order  to  download  the  plotting  data  file  from  the  VAX  to  an 
IBM  PC  the  user  must  first  logon  to  the  VAX  from  the  PC  using  the 
procedure  described  in  section  1.2  of  this  manual.  Once  the  user 
has  logged  on  to  the  VAX  he  should  then  switch  to  the  directory 
containing  the  sensitivity  algorithms  using  the  following 
command: 


$ SET  DEF  [.SENSITIV] 


The  user  should  then  switch  to  the  disk  drive  on  the  PC  where  he 
wishes  to  have  the  data  sent.  This  can  be  done  by  invoking  the 
DOS  shell  available  in  Procomm.  The  shell  is  invoked  by  simul- 
taneously depressing  the  Alt  and  F4  keys.  This  will  return  the 
user  to  DOS  and  allow  him  to  change  disk  drives.  For  example,  if 
the  user  started  Procomm  from  disk  drive  A and  wished  to  send  the 
plotting  data  to  an  existing  subdirectory  DATA  located  on  drive  C 
then  he  would  enter  the  following  sequence  of  commands: 


$ <Alt-F4> 


The  user  will  be  returned  to  DOS  and  see  the  following: 


ProComm  <==>  DOS  Gateway  Enter  EXIT  command  to  return 


The  IBM  Personal  Computer  DOS 

Version  3.10  (C)  Copyright  IBM  Corp  1981,  1985 

(C)  Copyright  Microsoft  Corp  1981,  1985 

A> 


The  user  should  then  change  disk  drives  and  directories,  and, 
return  to  the  VAX  using  by  entering  the  following  commands 
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A>C : 

C>CD\DATA 

OEXIT 


The  user  will  then  return  to  the  VAX,  having  changed  the  drive 
and  directory  on  the  PC. 

Once  this  has  been  accomplished  the  user  can  transfer  the  data 
from  the  VAX  to  the  PC.  To  begin  this  process  he  should  invoke 
KERMIT  file  transfer  protocol.  This  is  done  using  the  following 
command: 


$ KERMIT 


After  doing  so  the  following  information  should  appear  on  the 
screen. 


VMS  Kermit-32  version  3.0.052 

Default  terminal  for  transfers  is:  _TTP2 : 

Kermit-32> 


The  user  should  then  place  the  VAX  in  the  host  mode  with  PC 
control  by  using  the  SERVER  command. 


Kermit-32>SERVER 


After  entering  the  server  command  the  following  message  will 
appear. 


Kermit  Server  running  on  VAX/VMS  host.  Please  type  your 
escape  sequence  to  return  to  your  local  machine.  Shut  down 
the  server  by  typing  the  Kermit  BYE  command  on  your 
local  machine. 


The  user  can  then  initiate  file  transfer  by  simultaneously  pres- 
sing the  Alt-K  keys.  Procomm  will  display  a menu  giving  the  user 
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the  option  of  getting  or  sending  files  between  the  VAX  and  the 
PC.  To  download  a plotting  file  the  user  should  pick  option  1, 
GET.  Procomm  will  then  ask  the  user  for  the  name  of  the  file 
stored  on  the  VAX  which  he  wishes  to  download.  The  file  will  be 
stored  on  the  PC  under  the  same  name.  After  entering  the  file- 
name the  data  will  be  sent  from  the  VAX  to  the  PC.  After  the 
user  has  finished  downloading  files  he  should  press  Alt-K  and 
select  option  3,  FINISH.  This  will  exit  the  user  from  the  KERMIT 
file  transfer  protocol  and  return  the  user  to  the  '$'  prompt. 


5 . 2 Plotting  Sensitivity  Data 


Once  the  plotting  data  has  been  transferred  to  the  PC  the  user 
can  generate  plots  of  the  sensitivity  functions  using  the  BASIC 
program  PLOTSENS . BAS . 

To  begin  execution  of  the  plotting  program  the  user  should  place 
the  floppy  disk  containing  program  PLOTSENS. BAS  in  disk  drive  a 
and  type  the  following  command. 


A>BASICA  PLOTSENS 


This  will  start  execution  of  the  plotting  program.  The  plotting 
program  is  menu  driven  and  somewhat  self-explanatory.  The 
program  allows  the  user  to  plot  the  following  information: 


1.  Values  of  selected  system  variables  versus  time  (Time  domain 
sensitivity) . 

2.  Values  of  selected  system  variables  versus  front  wheel  steer 
angle  (Steady  State  sensitivity) . 

3.  Values  of  selected  elements  of  the  transfer  function  matrix 
versus  excitation  frequency  (Frequency  domain  sensitivity) . 

4.  First  order  standard,  percentage,  logarithmic  and  second 
order  standard  and  logarithmic  sensitivity  functions  versus 
time  (Time  domain) . 

5.  First  order  standard,  percentage,  logarithmic  and  second 
order  standard  and  logarithmic  sensitivity  functions  versus 
front  wheel  steer  angle  (Steady  state  sensitivity) . 

6.  First  order  standard,  percentage,  logarithmic  and  second 
order  standard  and  logarithmic  sensitivity  functions  versus 
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excitation  frequency  (Frequency  domain) . 


7.  General  sensitivity  measures  for  time  domain,  steady  state, 
and  frequency  domain  sensitivity. 


The  plotting  program  is  limited  to  ten  curves  per  graph. 

The  plotting  program  assumes  that  the  plotter  is  connected  to  the 
first  serial  port  (Coml) . If  your  plotter  is  connected  differen- 
tly it  will  be  necessary  to  edit  line  3970  of  the  program 
PLOTSENS . BAS  to  make  the  appropriate  changes. 

Once  the  user  has  quit  the  plotting  program  he  may  return  to  DOS 
by  entering  the  following  BASICA  command: 


SYSTEM 


5 . 3 Sample  of  Sensitivity  Plots 


Figures  2-7  represent  samples  of  graphs  obtained  for  the  2 
degree  of  freedom  nonlinear  spring-mass-damper  system  shown  in 
Figure  1.  This  system  has  been  widely  used  throughout  the  user's 
manual  to  explain  the  development  of  model  subroutines  for 
various  types  of  sensitivity  analysis. 

The  Figures  2-7  were  produced  using  the  plotting  program 
PLOTSENS . BAS . The  first  3 seconds  of  the  system  response  caused 
by  kinematic  excitation  (see  section  2.2.2)  is  shown  in  Figure  2. 
The  influence  of  various  system  parameters  were  determined  for 
the  following  system  variables:  the  vertical  positions  of  masses 
1 and  2 ( Z ^ and  Z2) , the  force  in  spring  and  the  force  in 
dashpot  . System  variables  and  Z2  represent  primary 
variables  (variables  obtained  directly  through  integration  of 
equations  of  motion)  while  forces  in  spring  K1  and  dashpot  C-, 
represent  secondary  variables. 

The  sensitivity  analysis  of  this  system  investigated  the  effects 
of  the  following  parameters:  the  springs'  stiffnesses  K-^  and  K-^ , 
the  dashpot  constants  and  ^(2,  and  the  damping  factors  Cn  and 
C2 . Parameters  K^,  K2 , |Hp  u 2 represent  primary  (constant) 
parameters  while  parameters  C-j_,  'and  C2  represent  secondary  (non- 
constant) parameters. 

Looking  on  Figures  2 - 7 it  can  be  seen  that  the  system  response, 
as  well  as  the  system's  sensitivity  functions  are  harmonic  in 
nature.  This  is  because  the  system  investigated  is  a purely 
vibrational  system. 
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The  system's  sensitivity  functions  show  the  effects  of  small 
parameter  changes  on  the  system's  response.  In  general,  the 
parameters  associated  with  sensitivity  functions  having  larger 
values  influence  a particular  variable  more  than  parameters 
associated  with  corresponding  smaller  sensitivity  function 
values.  This  allows  us  to  rank  parameters  in  their  order  of 
influence  on  system  response.  For  example,  Figure  3 shows  that 
at  t = 1.5  seconds  spring  stiffness  has  the  greatest  influence 
on  the  position  of  mass  m^.  At  t = 1.5  seconds  Figure  4 shows 
that  parameter  also  has  the  greatest  influence  on  the  position 
of  mass  m2 • Based  on  these  observations  we  can  conclude  that  the 
positions  of  masses  m^  and  m2  at  t = 1.5  seconds  can  be  effected 
more  by  changing  the  spring  stiffness  K-^  than  by  any  other 
parameter  considered. 

The  Figures  3 and  4 were  obtained  utilizing  the  percentage 
sensitivity  function,  which  is  one  of  several  types  of 
sensitivity  functions  which  can  be  used  and  are  described  in  this 
report.  Figure  6 represents  the  first  order  standard  sensitivity 
functions  for  all  6 system  parameters  considered  K-,  , K2 , C-^,  C2 , 
JU1  and  yu2 . Detailed  analysis  of  these  sensitivity  “results  should 
be  performed  in  conjunction  with  Figure  5 which  represents  the 
time  history  in  generating  the  spring  force  K-^,  due  to  frequent 
changes  in  spring  force  values  as  a result  of  system  vibrations. 

The  usefulness  of  the  general  sensitivity  measure  is  demonstrated 
in  Figure  7 in  which  the  percentage  sensitivity  function  values 
were  plotted  for  primary  variables  Z-±  and  Z2,  as  well  as  for 
secondary  variables  Fs  and  Fp.  The  general  sensitivity  measure 
indicates  the  global  effect  of  all  6 parameters'  changes  on 
system  variables.  For  example  as  seen  from  Figure  7,  during  the 
time  interval  1-2  seconds  the  order  of  influence  of  the  6 element 
parameter  vector  on  system  variables  is  as  follows:  spring  force, 
dashpot  force,  vertical  position  of  mass  m2  and  mass  m-^ 
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Figure  2. 


Vertical  Position  of  Mass  and  of  Spring-Mass- 
Dashpot  System. 
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Figure  3. 


Percentage  Sensitivity  of  Vertical  Position  of  Mass  m 
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Figure  5.  Forces  in  Spring  K1  and  Dashpot  Cl. 
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FIRST  ORDER  STANDARD  SENSITIVITY 
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Figure  6.  First  Order  Standard  Sensitivity  of  Spring  Force  Kl. 
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Figure  7.  General  Percentage  Sensitivity  Measure  of  System  Variables. 
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